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FOREWORD 

This  report  was  prepared  by  the  Metals  and  Ceramics  Division,  Materials 
Laboratory,  Air  Force  Wright  Aeronautical  Laboratories,  Air  Force  Systems 
Command,  under  Project  No.  2418,  ''Metallic  Materials,"  Task  No.  241803, 
"Behavior  of  Metals".  The  research  was  conducted  by  Dr.  David  I.  G.  Jones, 
Materials  Laboratory  (AFWAL/MLLN)  and  Dr.  hab.  Agnieszka  Muszyifska, 

Institute  of  Fundamental  Technological  Research,  Polish  Academy  of  Sciences, 
Warsaw,  Poland,  in  part  while  both  authors  were  at  Institut  National  des 
Sciences  Appliquees,  Lyon,  France,  during  the  period  November  1975  to 
August  1976.  The  subject  of  turbine  and  compressor  blade  vibration  pre¬ 
diction  and  control  is  of  great  importance  to  the  United  States  Air  Force, 
and  it  is  for  this  reason  that  a  complete  review  of  the  work  accomplished 
has  been  prepared.  The  report  was  typed  by  Judy  Mann  and  Kathy  Lee  Fox. 

This  report  covers  work  conducted  during  the  period  November  1975 
to  June  1979. 
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SECTION  I 
INTRODUCTION 

The  past  several  years  have  witnessed  renewed  interest  in  the  effects 
of  slip  on  dynamic  response  of  turbine  and  compressor  blades.  This  is 
partly  because  of  a  growing  recognition  that  although  energy  dissipation 
mechanisms  in  blade-disc  systems  are  still  not  well  understood,  this 
dissipation  is  one  of  the  most  important  factors  controlling  the  dynamic 
stresses,  and  partly  because  of  greatly  Improved  computational  capabilities 
which  have  permitted  linear  analysis  to  reach  a  quite  sophisticated  state 
of  development  and  hence  made  the  gap  in  our  knowledge  very  visible. 
Nonlinear  analysis  has  not  yet  reached  the  same  level  of  development 
(References  1-14). 

The  purpose  of  this  report  is  to  describe  a  combined  analytical/ 
epxerimental  investigation  of  a  compressor  blade  having  a  simple  dovetail 
root,  with  a  view  to  clarifying  the  effect  of  gross  slip  at  the  root  on 
the  dynamic  response  behavior  under  harmonic  excitation.  The  analytical 
part  of  the  investigation  examines  the  model ization  of  the  blade  in  terms 
of  a  simple  inertia-spring  system  with  coulomb  type  frictional  forces  at 
the  blade-disc  interface.  For  simplicity,  the  disc  is  assumed  to  be 
infinitely  rigid.  Of  course,  this  restriction  must  be  relaxed  before  a 
complete  blade-disc  assembly  is  analyzed,  but  we  are  not  at  this  point  yet. 

The  equations  of  motion  derived  on  the  basis  of  this  simple  model  are 
solved  by  a  method  of  harmonic  balance  assuming,  in  effect,  that  under 
cyclic  excitation  the  blade  will  exhibit  cyclic  response  at  the  same 
frequency.  The  solution  so  obtained  is  examined  numerically  to  determine 
regions  of  existence  and  nonexistence,  i.e.,  to  determine  the  frequency 
range  over  which  slip  occurs  and  that  over  which  it  does  not  occur 
(References  15-18).  The  agreement  between  analysis  and  experiment  is 
found  to  be  good  for  the  fundamental  mode  of  the  blade,  and  reasonably 
good  for  the  next  two  modes.  The  experiments  were  conducted  on  a  blade 
in  a  simple,  heavy,  fixture  with  the  normal  load  at  the  dovetail  provided 
by  two  thin  wires.  • 
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The  outcome  of  the  investigation  is  that  one  can  predict  with  a 
measure  of  accuracy  the  dynamic  response  behavior  of  a  simple  dovetail 
blade,  provided  that  the  appropriate  modal  information  is  either  measured 
or  predicted  by  more  detailed  analytical  procedures,  such  as  finite  element 
methods.  Extension  of  the  analytical  method  to  more  complex  blade-disc 
systems  is  a  task  for  the  future. 
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SECTION  II 

ANALYSIS  OF  RESPONSE  OF  A  BLADE  WITH  SLIP  AT  THE  DOVETAIL 

1.  MODELIZATIQN  OF  BLADE  IN  MODE  I 

The  type  of  blade  being  considered  in  the  analysis  Is  illustrated  in 
Figure  1.  We  shall  be  concentrating  attention  on  analysis  of  the  funda¬ 
mental  mode  vibration  response,  and  It  is  assumed  that  the  blade  resonant 
frequencies  are  well  separated.  The  blade,  when  vibrating  in  its 
fundamental  mode,  can  be  modelled  in  terms  of  the  simple  mass-spring 
system  shown  in  Figure  2,  where  the  inertias  1^  and  I2,  and  the  spring  K^* 
are  chosen  to  reproduce  the  observed  modal  stiffness  and  inertia  charac¬ 
teristics  of  each  individual  blade  geometry.  Ij,  I^,  and  K^*  will  vary 
slightly  with  the  rotation  speed  of  the  blade,  but  this  effect  is  neglected 
here.  It  can  be  allowed  for  in  the  case  of  any  real  blade,  but  would 
simply  entail  changing  the  parameters  I-j,  I^,  and  at  each  rotation 
speed  for  which  calculations  are  carried  out.  The  system  is  modelled  in 
terms  of  rotational  motion  because  the  type  of  motion  at  the  root,  when 
slip  is  occurring,  is  most  naturally  expressed  in  such  terms. 

2.  EQUATIONS  OF  MOTION 

For  the  nonlinear  case  where  slip  of  the  blade  root  takes  place, 
the  boundary  condition  seen  by  the  vibrating  blade  changes  from  an 
essentially  clamped-free  condition  to  an  essentially  pinned-free  condition. 
For  the  fundamental  mode  of  the  blade,  the  system  can  be  represented  with 
some  degree  of  accuracy  by  the  model  shown  In  Figure  2.  The  equations 
of  motion  of  the  system  with  slip  are: 

I101+K1  (1+ii^)  (©-,  -  02)  =  SLe1^1'  y)cos  a  (1) 

I292+Ki  (l+in1 )  (02  -  9})  +  (uNR0/cost/>)sgn  e2  =  0  (2) 


3 


AFWAL-TR-80-4003 


The  term  wNRQ/cos^  is  the  frictional  moment  opposing  the  motion.  Figure  2 
shows  the  geometrical  factors  involved  in  this  derivation,  a  is  the  angle 
between  the  blade  direction  at  the  tip  leading  edge  (where  S  is  applied) 
and  the  root  axis  as  illustrated  in  Figure  1.  From  vertical  equilibrium 
in  Figure  2: 

2F^cosi|i  =  N  i.e.  F-j  =  N/2cost)i 
F2  =  ^1  3  PN/2cost)< 

Moment  =  ZF^Rq  =  pNRQ/cos^  (3) 

Equations  1  and  2  can  be  rewritten  in  the  form: 

I1(e1/S)+K1(l+in1.)  (e^S-Gj/S)  =  e1' Ut~Y)Lcosa 

I2(02/S)+K1(l+in1)  (02/S-91/S)+(uNRo/Scos^)sgn  02  =  0 

This  shows  that  the  mobility  0^/S  depends  on  pN/S,  so  that  the  analysis 
within  its  limits  is  valid  for  low  or  high  values  of  S  or  of  pN;  it  is 
the  ratio  of  these  two  quantities  which  is  important. 

3.  SOLUTIONS  FOR  STEADY  STATE  MOTION 

a.  Solution  for  Non-Slip  Condition 

If  S/pN  is  sufficiently  small,  slip  cannot  occur  and  Equations  1 
and  2  reduce  to: 


I1 91+K1  ( l+1n] )  (Oj)  =  Se'l^t“y^Lcosa 
the  solution  of  which  is: 

X1  =  0 1 L  =  SL2cosael(ajt'Y)/[-I1u)2+K1(l+in1)] 


(4) 


(5) 


5 


AFWAL-TR-80-4003 


The  values  of  I-j,  ,  and  n-j  needed  to  model  a  specific  system  are  now 
determined  by  comparing  the  response  X-j  with  the  measured  response  for  low 
values  of  S/yN.  The  measured  resonant  frequency  gives  one  such 
relationship 

I^2  *  K1  (6) 


The  measured  response  at  resonance  gives  another: 

lx 


'1 


2  2 

tu,  L  COSci  ,  2„_ 

1 _  _  L  cosot 


Kin, 


I 


lnl 


:.  I,  =  Sl/cosct/r^  |x] 


(7) 


and  the  measured  damping  n-j  is  given  by: 

Auj-j 


Where  Au>-|  is  the  difference  between  the  two  frequencies  at  which  the 
response  is  3  db  below  the  peak  response  (3.01  actually.'). 

For  example,  in  the  case  of  the  blade  used  in  the  experimental 
investigations,  such  calculations  are  summarized  in  Section  V  and  the 
average  parameters  were  determined  to  be: 

^  =  1.04  x  10'3  kg  m2 
=  635  Nm/radian 
n1  =  0.010 
w1  =  124.2  Hz 

for  the  fundamental  mode. 
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b.  Solution  with  Slip 

For  simplicity,  we  let  n-|  =  0  in  Equations  1  and  2;  for  the 
purposes  of  this  part  of  the  analysis  it  is  sufficient.  If  we  seek  a 
steady  state  solution  of  these  equations,  then  to  the  first  order  we 
assume  that: 

02  =  ©2  sin  cot 

©1  =  9-|-]  sin  cot  +  ©i2  cos  cot  (9) 

S  =  S-j  sin  cot  +  S2  cos  cot 

Furthermore,  we  consider  only  the  first  term  in  the  expansion  of  sgn  (02) 
{=  sgn  (cos  cot)}  in  ascending  terms  of  sin(ncot)  and  cos  (not),  so  that: 

sgn  (cos  cot)  ?  ^  cos  cot  (10) 

Then,  introducing  these  equations  into  Equations  1  and  2,  and  comparing 
terms  in  sin  cot  and  cos  cot,  we  have  the  four  algebraic  equations  for 

^11  *  *12  ’  ~*1  *  ^od  $2  ■ 

( - 1  ^  co^  )®ii“K-|©2  =  S-|Lcosa 
(K^-I-|Co^)0i2  =  S2Lcosa  (11) 

(-I2co2+K1)02-K1011  =  0 
-K-|0-|2  =  4pNR0/ircos^ 


The  solution  of  this  equation  is  readily  obtained.  From  the  values  of 
0^ ,  0^2>  S.| ,  and  S2  so  determined,  we  then  can  write  the  solution  for 
0^  and  S  in  this  form: 


61  = 


-S-j  (K-j-^co  )Lcosasin  cot  / 


4pNR 


\  cOSijmK 


t) 


cos  cot 


(12) 


2  \  /■.  tin 


SLcosa  =  S^Lcosasincot-(K^-I^co  )  (4pNR0/TtK^costJc)coscot 


03) 
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Clearly,  we  can  eliminate  in  these  two  equations.  This  is  most 
simply  done  for  our  present  purposes  by  first  calculating  J 9 ^  |  and  |S|, 
since  these  are  the  quantities  we  would  ordinarily  measure.  Then: 


S12(K1-I2oj2)2L2cos2ct  /  4pNRQ 
u>4  +  2tJj2 ^2  C0S'l'1Tl<'i 


(14) 


L2cos2ot | S |  2  =  S^2Lcos2a+(Ki-Iito2)2(4uNRo/7rKiCosiJ;)2 


(15) 


Eliminating  S-|  and  simplifying  gives: 


+ 


2  4  2 

B-j  L  cos  a 


(16) 


L4cos2a[l  (1-I1tu2/K1  )2] 
K12(I1oJ2/K1)2[l+I2/I1-(I2/I1)(I1aJ2/K1)  ]2 


(17) 


where 


8^  =  4uNRQ/TrSLcos^cos  a  =  ( 4/tt )  S2 


(18) 


B2  =  uNRo/SLcos^cosa 


If  8-j  =  0,  as  for  a  freely  pinned  blade  or  for  S  very  large,  then  |X^| 
has  an  infinite  amplitude  at  the  frequency  where: 


Ip  I i 

(i  +r)r  =  r+ 1 

M  l2  l2 


(19) 
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This  allows  us  to  estimate  I2,  since  the  pinned-free  mode  of  the  particular 
blade  used  was  determined  experimentally  to  be  429  Hz. 

-p-  +  i  =  (42-9)-  J  =  11.93 
l2  (124.2)^ 

I2  =  9.51  x  10'5  kg  m2 

for  the  particular  blade  used  in  the  experimental  investigation.  A  cal¬ 
culator  program  to  determine  J X^/S  |  and  |X^/S|  is  given  in  Table  1.  It 
is  seen  from  Equation  17  that  |X2/S|  can  exist  only  over  a  narrow  range 
of  frequencies,  and  it  is  only  within  this  band  that  slip  can  occur. 

Outside  this  range,  although  |X,/S|  as  given  by  Equation  16  may  exist, 
the  original  assumption  that  |X2/S|  exists  is  violated  and  hence  we  must 
return  to  the  linear  solution  given  by  Equation  5.  Figures  3  to  10  show 
the  calculated  values  of  |X^/S|  and  (X2/S|  versus  frequency  for  B2  =  50, 

20,  10,  5,2,1,tt/4,  0.3  and  for  the  aforementioned  values  of  1^,  I2>  , 

and  .  Figure  11  shows  a  summary  of  the  solution  for  several  values 
of  02<  These  analytical  results  will  be  compared  with  experiment  in 
Section  V. 

c.  Graphical  Solution  for  Slip  Motion  and  Hysteretic  Damping 

The  equation  of  motion  of  the  system  with  slip  and  hysteretic 
stiffness  can  be  written  in  the  more  general  form: 

••  k-iDl  •  1  • 

I101  +  -^-L  9-,+^e.j-  j  (-^  92  +  k102)(l+sgn  3)  =  S^osuit  (20) 

_  •  k-.ru  '  '  1 

[I202  +  yNR  sgn  «2-  (e1-e2)-k1e1]  2  (1+sgn  8)  +  k102  =  0 

where 

S^  =  S  cosa 


9 


AFWAL-TR-80-4003 


TABLE  1 

HP67  CALCULATOR  PROGRAM  TO  CALCULATE  RESPONSE 

Input;  put  in  in  x  register  (Hz) 


Manual  storage: 


Storage 

II 


II 


location  0 
"  1 

"  2 

"  3 

"  4 

"  5 

"  6 


Kx (635  Nm/rad) 
ui  (120  Hz) 

1 2/1 1  (.0915) 

8 2  UO) 

n !  (.01) 

L ( . 20) 
cos  a  ( . 724 ) 


Output:  Label  A  -  ] X j /S |  linear  (Equation  5) 

Label  B  -  | Xj/S 1  slip  (Equation  16) 

Label  C  -  |x2/S|  slip  (Equation  17) 

Note:  Label  A,B,  and  C  must  be  operated  in  sequence. 


Internal  storage 

Storage  location  A  -  w2  rad2/sec2 


II 
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Figure  3.  Predicted  Response  of 

Blade  in  Mode  1  for  62  =  50 
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Figure  4. 
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Figure  5.  Predicted  Response  of 

Blade  in  Mode  1  for  =  ^ 


Figure  6.  Predicted  Response  of 

Blade  in  Mode  1  for  B2  =  J 
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FREQUENCY-  Hz 


Figure  11.  Graphs  of  ]X^/S|  and  | Xg/S |  versus 
Frequency  (Summary) 
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For  8  <  0  :  02  =  0  -  one  degree  of  freedom 
For  8  >  0  :  two  degrees  of  freedom 

8  =  0  is  tne  condition  of  the  threshold  of  slip  (it  will  be  defined  later) 
R  -  geometrical  parameter  (RQ/cos^) 
n-j-  hysteretic  damping  coefficient 
The  solution  of  Equation  20  is: 

02(t)  =  A  cos(cot  +  y)  (21) 

put  into  Equation  1  and  expanding  the  function  sign  [-Au)Sin(u)t  +  y)l 
in  the  Fourier  series  and  leaving  only  the  first  term  we  reduce  the 
problem  to  solving  the  linearized  system 

1"  k,rN 

I101  +  —  (0t_92^  +  (0i'02^  =  S]C0SCl)t 

i2e2 +  -ir  ^2~ei^ +  ki  ^02_0i^ =  S9nA  (22) 

looking  for  the  solution  in  the  form 


0-j (t)  =  D  cos  (ut  +  y-6) 
0„(t)  =  A  cos  (wt  +  y) 

i  <- 


(23) 


if  the  condition  8  >  0  is  satisfied.  If  we  have  8  £  0  we  look  for  the 
solution 

0^(t)  =  D*  cos  (wt-Y*)  (24) 


of  the  equation 


k,r)-i 

1,0,  +  8,  +  k,0,  =  S,cos  ut 

11  u  1  II  I 


(25) 
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The  last  one  gives 


D* 


^1 

,  6*  =  arctan — , 

k1  -^ur 


(26) 


Putting  Equation  23  into  Equation  22  we  get 

A^u%.|  -  2Au)2c^  (sgn  A ) +  ^[(r^v)2  +  (v-1)2]  =  S2  (27) 

❖l  =  [I2(v-1  )-Ij  ]  +  (n-j vl2) 

<t>  2  =  n1v[2I2(v-l)-I1] 

4>3  =  {(v-l)[I2(v-l)-I1]  -  I2(niv)2}2 

I^a)2 

v  =  - o~ 

k1(l+n1  ) 
a,  =MR 

I  IT 


From  Equation  27  we  set  the  amplitude  A: 


s, 

|  /5~3  Si 

sl  1  /4>  3 

5T<  h 7 

—  (v-1)  z  +  (niv)  z 

«1 

nr~^7~/{  1)Z  +  (niv) 

A= 

to  ^ 

2  "”2 

n- - z — 

±  +  /»!!  i~a  1  <f>  3 

<*>2$  1 

p - 2 — 

A=±a  i<f>2  —  ^^l^l- "0ll4,3 

“)zif  1 

amplitude  does 
not  exist 

where : 

“2  -  k'  (1+,'i|  [2T7  +  I,] 

and  the  phase  angle 


Awzn 1 I2 v  +  (l-v)aj 

y  =  arctan  - — - 

aim v-Aw2  [I2(v-1)-Ii] 
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The  amplitude  of  the  solution  (Equation  23)  is  given  in  the  table  on  the 
previous  page.  As  one  can  see  for  some  values  of  the  parameters  we  have 
two  values  of  A  (one  positive  and  one  negative)  or  four  values  (real  ones!) 
In  some  region  the  solution  (Equation  23)  does  not  exist;  it  is  the 
domain  of  the  solution  (Equation  24)  of  the  reduced  system. 

The  coefficient  of  the  threshold  of  slip  "6"  can  be  defined  as  shown 
in  the  above  table: 

si  1/T3  "Sa  |  (v-1)  [I2  (v-D-Iil-Ii  (niv)2| 

o  1/  ' 


S1  ,/T?  *5a  J  (v-1) [I2 (v-l)-I 

a!  »  4>  j  4uNR  /[  I2  (v-1)  -Ii)2  + 


The  relationship  between  A  and  the  amplitude  D  is 


(n,vl2  )  2 


7+  /..(■-  ¥)'•  hs; 

(29) 

/l+m  V 

An  j i2u2  +  0] 

arctan  - 

2 

(30) 

AI 2 w 2 -Ak j  (l+m)-niai 

where  A  is  calculated  from  Equation  27.  (See  table) 

A(oj)  curves  are  plotted  graphically  as  in  Figure  12 


k i (1+n i ) 

'  1 

2 

]F7r 

r 

2 (1-n  i ) 

.  I2 

II 

'i22  ix 

1 2 / 

l+m  1  Ij 


Different  levels  of  — give  different  shapes  of  the  |A|(u>)  curve 
al 

(qualitative  shape!) 
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Figure  12.  Graphs  of  |Aj  and  S./a,  versus 
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The  case  n-|  =  0  (Figure  13)  gives: 


(31) 


A  method  of  connection  of  the  solution  D  from  A  is  given  in  Figure  14. 

The  curves  are  plotted  in  an  acceleration  scale  for  comparison  with 
the  experimental  data.  It  is  only  a  qualitative  picture  here. 

4.  CALCULATION  OF  APPARENT  DAMPING 

As  Figures  3  to  10  show,  the  response  mobility  | X-j /S |  of  the  blade 
does  not  have  the  classical  damped  single  degree  of  freedom  system  shape 
but  is  sharply  "cut-off"  when  slip  occurs.  Strictly  speaking  the  system 
is  not  "classically"  damped  and  cannot  be  assigned  a  true  measure  of 
modal  damping.  However  the  shape  of  |X^/S|  as  a  function  of  frequency 
depends  only  on  and  for  each  value  of  3^  this  shape  is  unique.  One 
may  therefore  define  an  "apparent  modal  damping"  by  the  usual  "half  power 
bandwidth"  definition,  i.e.,  in  Figures  3  to  10  we  seek  the  frequencies 
fj  and  f^,  when  the  response  is  1//?  times  the  peak  amplitude  (-3.01  db) 
and  define  the  apparent  damping  as: 


9 


(32) 


On  this  basis,  we  can  plot  a  graph  of  versus  1 / 62  as  in  Figure  15. 
H-.  increases  as  l/Bo  increases,  at  least  within  the  range  of  values 

^  t 

examined  here.  For  l^lO.Ol,  the  value  of  n-j  remains  at  the  non-slip 
value  n-j  -  0.01  observed  in  the  tests. 
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Figure  14.  Graphical  Method  of  Connecting  |D|  and  |A! 
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One  can  also  examine  the  bending  moment  at  the  root  as  a  function 
of  &2-  The  value  of  M  is  given  by: 


M  = 


K,e, 


M 

K] 

*1 

SL 

L2 

S 

(33) 


For  a  specific  blade,  one  can  determine  the  relationship  between  M  and 
the  maximum  stress  in  the  fundamental  (or  any  other)  mode,  but  at  present 
it  is  sufficient  to  examine  | M/SL |  as  a  function  of  3 ^  as  determined  from 
Figures  3  to  10.  The  resulting  graph  of  | M/SL |  versus  ^  shown  in 
Figure  16. 


5.  CALCULATION  OF  EFFECTS  OF  ROTATION  SPEED  ON  APPARENT  DAMPING 

It  is  generally  difficult  to  predict  the  response  of  a  blade  at 
high  rotational  speed  because  the  exciting  forces  on  a  blade  are  not 
known  a-priori.  However,  for  a  given  blade  geometry  S  will  depend  on  the 
rotation  speed  Q  in  the  manner: 


(34) 


where  Sf)  is  the  force  at  some  speed  :iQ  and  n  is  an  exponent  not  neces¬ 
sarily  equal  to  2.  S  will  be  a  maximum  at  the  various  multiples  of  the 
blade  passage  frequencies,  so  that  the  excitation  force  at  any  condition 
where  a  blade  passage  frequency  is  near  to  a  natural  frequency  will  be 
nearly  harmonic. 


The  normal  force  N  resulting  from  the  rotation  speed  .<  is  given, 
as  a  crude  approximation,  by  a  relationship  of  the  type: 


I'TS/Wj 
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Figure  16.  Graph  of  Predicted  Value  of  | M/SL |  versus  &2  for  Mocle  1 
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where  M  is  the  blade  mass,  R  an  average  radius  and  :l  the  rotation  speed 
in  rpm.  It  is  seen  that  if  n  =  2,  &2  is  constant  at  all  speeds  and  hence 
that  the  damping  due  to  slip  will  be  the  same  at  all  speeds,  which  is 
contrary  to  observation.  Hence  one  must  conclude  that  in  most  cases 
n  <  2.  For  example,  if  we  assume  that  n  =  1.5: 

1.5 


3, 


?  n  5 
4tt  pMRR  ft 
o  o 


2  2 
SoLcosi^cosa(60) 


(36) 


As  a  hypothetical  example,  consider  the  cases  where  u  =  0.15,  0.1  and 
.,.05,  M  =  0.235  kg,  R  =  0.5m,  Rq  =  0.01m,  =  15000  rpm,  SQ  =  40N, 

_  =  0.2m,  cosij,  =  0.5,  cosa  =  0.724.  Then: 

B2  =  8. Op  Q°‘ 5  (37) 

For  each  value  of  p  and  U,  one  can  calculate  B2  and  hence  read  the  value 
of  effective  loss  factor  off  Figure  15.  Some  predicted  graphs  of  the 
effective  loss  factor  versus  Q  are  plotted  in  Figure  17.  It  is  seen 
that  for  p  =  0.15,  the  behavior  is  very  similar  to  that  observed  in  a 
test  engine  staoe  (Reference  2),  but  that  as  y  is  decreased,  the  apparent 
modal  damping  increases  rapidly.  It  is  seen  from  Equation  36  that  it  is 
desirable  to  reduce  p,M,R  and  RQ,  and  increase  So>  L ,  cost*;  and  cosa,  as 
far  as  possible,  to  achieve  high  "damping". 

6.  M0DELIZAT 10N  OF  BLADE  IN  SECOND  BENDING  MODE 

For  the  second  bending  mode  of  a  blade,  or  for  any  mode  in  which 
the  deformations  are  primarily  bending,  the  model  selected  for  the  funda¬ 
mental  mode,  illustrated  m  Figure  2,  is  directly  applicable  and  only  the 
value  of  1^,  I9  and  K^*,  will  be  different.  This  means  that  Equations  1 
and  2  and  the  solutions  (Equations  16  and  17)  are  directly  usable.  From 
the  low  level  tests  described  in  Section  V,  the  corresponding  values  of 
I, ,  K.  ,  and  are: 

iii  -A  2 

^  =  1.62  x  10  Jkg  m 
K]  =  11000  Nm/rad 
n-j  =  0.0067 
oj-|  =  415  Hz 
I2/I1  =  0.150 


24 


AFWAL-TR-80-4003 


Figure  IF.  Typical  Predicted  Variation  of  with  Q 

The  corrt  ::'"Miing  solution  for  ^  defined  as  before,  equal  to  50,  20, 
10,  5,  2,  1,  74,  and  0.3  are  illustrated  in  Figures  18  to  25. 

The  difference  between  the  second  mode  solution  and  the  fundamental 
mode  solution  is  seen  to  result  only  from  the  altered  values  of  I-j,  K1  , 
n1  ,  and  .  This  will  be  true  for  all  the  essentially  "bending-type" 
modes  for  which  the  moment  at  the  root  caused  by  the  blade  deformation 
is  in  one  direction  only,  and  is  not  a  twisting  type  of  moment,  as  for 
the  torsional  modes. 

The  value  of  is  determined  from  the  ratio  of  the  second  mode 
resonant  frequencies  in  the  clamped-free  and  pinned-free  conditions. 

For  this  blade  =  2.43  x  10  ^  kg  m^. 
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Figure  18.  Predicted  Response  of 
Blade  in  Mode  2  for  B9 


Figure  19. 
50 


Predicted  Response  of 
Blade  in  Mode  2  for  ~ 
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Figure  20.  Predicted  Response  of 
Blade  in  Mode  2  for  0^ 


Figure  21.  Predicted  Response  of 
10  Blade  in  Mode  2  for 
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Figure  22.  Predicted  Response  of 

Blade  in  Mode  2  for  82  =  2 


Figure  23.  Predicted  Response  of 

Blade  in  Mode  2  for  62  = 
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Figure  24.  Predicted  Response  of 
Blade  in  Mode  2  for  B9 


Figure  25.  Predicted  Response  of 
tt/4  Blade  in  Mode  2  for  8,  = 
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7.  MODELI ZATION  OF  BLADE  IN  FIRST  TORSIONAL  MODE 

The  third  mode  observed  in  the  response  spectrum  appears  to  be  a 
torsional  type.  For  this  case,  the  moment  at  the  root  when  slip  occurs 
is  of  a  twisting  type,  as  illustrated  in  Figure  26.  In  this  case  I-j 
oecomes  the  torsional  inertia  of  the  blade,  I2  that  of  the  root,  and  K-j 
the  torsional  stiffness.  The  equation  of  the  blade  therefore  becomes,  in 
any  predominantly  torsional  type  of  mode: 

1 1  tj-j +K.|  (1  +  ir^J  =  Sbe^wt’Y^cosa 

1 2®2+K7  (1+iri1)  (O2-01  )  +  (uNbo)sgn02  =  0 

when  0-|  and  02  are  now  the  torsional  deflection  of  the  blade,  i.e.,  for 
this  type  of  mode  X-j  =  0-jb  and  X-,  =  02b. 

I,  (0,/S)+K1  (l+ir)l)(0,/S-09/S)  =  bcosa  ei(ajt-^ 

11  1  112  (38) 

I2(02/S)+K1  (1+in, )(02/S-01/S)  =  (uNb0/S)sgne2 

We  tnerefore  recover  the  fundamental  mode  solution  of  Equations  16  and  17 
if  we  replace  L  by  b  and  uNRQ/SLcosa  cosi{>  by  uNbQ/Sbcosa.  The  correspond¬ 
ing  values  of  I-j,  I2,  ,  and  w.  are  given  in  Section  5.  Figures  27 

to  34  show  some  typical  calculated  results. 

8.  NORMAL  MODE  ANALYSIS  OF  BLADE  RESPONSE  WITHOUT  SLIP 

In  order  to  establish  the  low  amplitude  level  behavior  of  the  blade, 
and  to  identify  parameters  in  a  preliminary  fashion,  it  is  essential  to 
first  consider  the  behavior  of  the  blade  with  no  slip  at  the  root.  In 
this  case,  in  the  fundamental  mode,  the  blade  behaves  essentially  as  a 
clamped-free  beam.  Since  the  blade  is  twisted,  it  is  necessary  to  define 
a  coordinate  system  within  the  middle  surface  as  illustrated  in  Figure  35. 
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Figure  29.  Predicted  Response  of 

Blade  in  Mode  3  for  Bo  =  10 


Figure  30.  Predicted  Response  of 
Blade  in  Mode  3  for  3 
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If  the  blade  is  then  subjected  to  an  exciting  force  S(t)  =  S6(x-x^)  6 
(y-y-|)sin  ut  normal  to  the  middle  surface  at  the  point  (x^y^),  then  the 
equation  of  motion  of  the  system  can  be  written; 

L3A(K*W)  +  u(32w/3t2)  =  Sfifx-x^  5  (y-y-^sin  ut  (39) 

where  K*  is  the  complex  stiffness  of  the  blade,  assuming  hysteretic 
damping  and  A  is  a  differential  operator  in  the  space  variables  x  and  y, 
and  u  is  the  mass  per  unit  area  at  the  point  x,y.  Since  the  blade  is 
usually  very  stiff  in-plane  and  much  more  flexible  normal  to  the  middle 
surface,  W  is  practically  normal  to  the  middle  surface  quite  frequently. 
Now  if  the  normal  modes  are  known,  from  experiment  or  analysis,  the 
displacement  W  can  be  expressed  as  a  series  of  normal  modes: 

00 

W=  I  W  <p  (x,y)  (40) 

n  =  1  n  n 

where  Wn  is  a  scalar,  depending  on  time,  representing  the  amplitude  of 
the  nth  normal  mode  and  <|>n  is  a  vector  representing  the  normalized  dis¬ 
placement  function  in  this  mode.  By  definition,  the  normal  modes  are 
the  set  of  discrete  functions  which  satisfy  the  homogeneous  equation: 

A(K$n)  -  (p  u2/L3Hn  =  0  (41) 

This  classical  relationship  allows  one  to  convert  Equation  39  from  a 
complex  partial  differential  equation  in  the  space  variables  to  an 
infinite  series  (hopefully  convergent)  of  algebraic  equations:  Putting 
Equations  40  and  41  into  Equation  39  gives;  in  complex  notation: 

00 

X  [(l+in)pun2Wn  +  yWn]d>n  =  Sexp(iut)6(x-x1 ) 6(y-y} )  (42) 

n  =  1 

one  must  now  make  use  of  the  orthogonal  property  of  the  normal  modes: 

/  A4ntmdxdy  =  6mn  (43) 

x  y 
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where  6 _ =  0  for  m  i  n  and  6^  =  1  for  m  =  n.  Therefore,  the 

mn  mm  * 

factoring  Equation  42  by  and  integrating  over  the  surface  gives: 

[ ( 1  +in)uj^W1+W1  ^u^.^dxdy  =  £  p. ^(x-Xj  )6(y-y1  )eiajt  dxdy 

where  <J>,  .<J>i  is  the  scalar  product  of  the  vector  4>^  with  itself.  This 
equation  can  clearly  be  written: 

mj[<^(1+io)-w^]W.|  =  f^Se^t  (44) 

t 

where  is  the  mass  of  the  blade  and  is  a  nondimensional  parameter: 

t 

m-j  =  /  /pdxdy  (45) 

x  y 


r 


/  /  n.|16(x-x1 )6(y-y1 )dxdy 


/  /  p<£,  .<j>.dxdy 
x  y 


m1n(x1  ,y1 ) . d>1  (x]  ,y] ) 

/  /p^.cLdxdy 
x  y 


(46) 


where  n  is  the  unit  vector  normal  to  the  middle  surface  of  the  blade  and 

t  t 

S  =  Sn(x,y).  If  we  let  m^  =  m^/8^,  this  means  that  the  behavior  of  the 
blade  in  the  fundamental  mode  is  equivalent  to  that  of  a  single  degree  of 
freedom  system  of  mass  ,  stiffness  =  m^co^,  and  loss  factor  n-  Note 
that  B-|  depends  on  the  point  x-j  ,y^ .  We  can  define  a  parameter  B-j  by 
letting  (x^  ,y^ )  =  (0,0) ,  the  tip  of  the  blade  at  the  leading  edge.  Then: 

B1  =  61  n”( x^  ,y] )  (x^  ,y] )  (47) 
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provided  that  |^(o,o)|  =  1,  i.e.,  the  modal  function  is  normalized. 
The  solution  of  Equation  44  then  becomes: 

S01  (x,y)eiwt[n(x1  ,y1 )  .4>1  (x^)] 

W  =  - - — 2 - ; - 2 - ' —  (48) 

m^ (l+in)-uj  ] 

For  excitation  at  (o,o)  and  measurement  of  the  response  at  (x,y),  we 
then  have  at  resonance: 

S|  <J>,  (x,y)  | 

|W|  =  - 1  2  -  (49) 

I  inm-|Uj.|  | 

•*k1(x,y)|  -  (|||n)  (50) 

This  means  that  we  can  measure  the  normal  mode  (or  any  other  mode) 
by  determining  | W/S |  for  excitation  at  any  fixed  point  and  pickup  at 
the  point  x,y  and  then  multiplying  by  n*  If  n  is  not  constant  for  all 
measurements  (as  a  result  of  extraneous  sources  of  hysteretic  damping 
such  as  cables),  then  Equation  50  allows  one  to  compensate  for  this 
source  of  error.  Note  that  n  can  be  measured  by  the  "half-power 
bandwidth"  method,  so  that  one  can  actually  determine  m^  directly  from 
the  peak  value  of  | W/S |  at  resonance,  using  Equation  49,  if  the  measure¬ 
ment  chain  is  properly  calibrated. 
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SECTION  III 

ANALYSIS  OF  BLADE  RESPONSE  WITH  SLIP  AT  A  PLATFORM 
1 .  INTRODUCTORY  REMARKS 

Compressor  and  turbine  blade  failures  caused  by  excessive  vibration 
can  and  often  do  arise  in  jet  engines  whenever  high  flow-induced 
exc’tation  forces,  high  static  stresses,  and  low  modal  damping  levels 
occur  at  the  same  time.  While  blade/disk  interactions  will  complicate 
the  response  behavior  of  eacn  individual  blade,  and  cause  circumferentially 
changing  peak  stress  levels  around  the  disk,  it  is  slip  at  the  blade/disk 
interface  which  provides  a  major  mechanical  source  of  damping,  in  addition 
to  the  aerodynamic  and  material  sources.  Attempts  to  increase  slip  damp¬ 
ing,  bv  means  of  mid-span  or  tip  shrouds,  or  by  means  of  mechanical  con¬ 
nections  between  adjacent  blades,  have  not  been  very  successful,  perhaps 
because  of  the  tendency  for  corresponding  points  on  neighboring  blades  to 
vibrate  with  only  relatively  small  amplitude  and  phase  differences  except 
where  extreme  efforts  are  made  to  mistune  adjacent  blades  relative  to  each 
other,  i.e.,  the  effect  is  usually  to  stiffen  rather  than  dissipate  energy. 
Certainly,  the  analytical  difficulties  of  predicting  the  response  of 
complete  blade/disk  systems  with  slip  at  each  blade/disk  interface,  or 
between  each  blade,  are  formidable  (References  12,  19)  and  will  not  be 
addressed.  We  shall  examine  a  configuration  in  which  each  blade  will 
slip  relative  to  the  disk  rather  than  relative  to  neighboring  blades. 

With  proper  attention  paid  to  the  static  and  dynamic  forces  involved,  such 
a  configuration  can  lead  to  high  slip  damping  even  at  high  rotational 
speeds.  The  response  behavior  of  such  a  system  is  highly  nonlinear,  so 
we  shall  assume  that  the  compliance  of  the  disk  is  infinite,  as  a  first 
approximation,  in  order  to  make  the  analysis  more  tractable.  Figure  35 
shows  such  a  blade  concept,  as  compared  with  a  simple  dovetail  root  and 
a  Christmas  tree  root.  The  gap  between  the  outer  step  of  the  blade  root 
and  the  disk  is  very  important,  since  it  must  close  only  at  the  selected 
rotational  speed  above  which  some  damping  is  reguired. 
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2.  ANALYSIS  Of  DYNAMIC  RESPONSE 

In  order  to  model  the  dynamic  response  behavior  of  the  blade  in  a 
single  mode,  usually  the  first,  different  physical  models  are  required 
depending  on  whether  slip  is  or  is  not  occurring  at  the  root  or  at  the 
sub-platform.  The  blade  is  represented  for  purposes  of  analysis  as  shown 
in  Figure  36.  The  various  physical  models  which  represent  the  blade  under 
different  conditions  are  shown  in  Figure  37.  The  parameter  6  defines 
whether  slip  occurs  at  the  lower  dovetail  or  whether  it  is  locked  at  this 
point,  whence  =  0;  and  6-|  defines  whether  the  mass  m^,  representing  the 
sub-platform,  is  slipping  against  the  disk  or  is  locked  in  place  by  the 
frictional  forces,  for  which  case  X^  =  0.  The  masses  m^ ,  m^,  and 
and  the  stiffnesses  k^  and  k^  must  be  selected  in  accordance  with  the 
blade  geometry.  The  disk  impedance  is  assumed  to  be  much  greater  than 
that  of  the  blade,  for  simplicity. 

By  experimental  or  analytical  (e.g.  finite  element)  methods,  one 
can  determine  the  ratio  of  the  response  at  any  point  j  to  the  force 
applied  at  any  point  i,  i.e.,  the  compliance  G^lcu),  and  this  data  is 
then  used  either  directly  in  a  modal  analysis  or  indirectly  to  determine 
m^  ,  m^,  m-j,  and  k^  ,  k^  for  the  discrete  element  model.  Both  methods 
should  give  comparable  results  if  the  respective  assumptions  and 
simplifications  are  consistent,  but  the  discrete  element  model  is  the 
easiest  to  analyze.  For  the  model,  in  Figure  38  the  equation  of  motion 
can  be  represented  for  all  cases  by  the  equations: 


=  S  cos  »t 


37 


AFWAL-TR-80-4003 


|m2X 


*]| 


X2  +  uNR  sgn  X2  +  k  (X2  -  Xi)  -  X  i  j  (1  +  sgn  6). 

sgn  (fi  -  J20)^  +  4k  X2  =  0  (52) 


I1' 


jm3X3  +  pNiRi  sgn  X3  +  k^  (X3  -  X3)  -  xj  +  k2^  X3j,  (1  +  sgn  Bj! 


+  2  (kj  +  k2)  X; 


1  +  sgn  (Q  -  S2  o ) 


=  0 


where 


k  =  kjk2/(k2  +  k2) 


(53) 


(54) 


is  the  equivalent  stiffness,  m, ,  m^,  and  are  masses,  k-j  and  are 
sub-stiffnesses,  q  is  the  blade  loss  factor,  S  and  to  are  the  amplitude 
and  frequency  respectively  of  the  exciting  force,  p  is  the  dry  friction 
coefficient,  N  and  are  the  normal  loads  between  the  blade  root  and 
the  disk,  R  and  R-j  are  coefficients  depending  on  the  root  geometry 
(References  17,  18)  8,  8-|  are  coefficients  of  the  slip  thresholds.  If 

we  let  Q.  be  the  rotational  speed  at  which  the  gap  A  -*•  0,  then  different 
solutions  occur,  depending  on  whether  Q<Qq  or 


a.  Case  1 ,  51  <  ilQ,  8  >  0 

The  analysis  is  the  same,  essentially,  as  that  for  a  blade  with 
a  simple  dovetail  root  (Equations  17,  18). 

Equations  51  to  53  reduce  to: 

mixi  +  —  (Xx  -  X2)  +  k  (X 2  -  X2)  =  S  cos  lot  (55) 


m2X2  +  y NR  sgn  X2  +  (X2  -  Xj)  -  kXz  +  2kX2  =  0  (56) 

We  look  for  an  approximate  steady  state  harmonic  solution  of  the  nonlinear 
Equations  55  and  56  in  the  form: 


Xx  (t)  =  D  cos  Ut  +  y  -  6)  and  X2  (t)  =  A  cos  ( cot  +  y)  (57) 
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This  is  accomplished  by  expanding  the  function  sgn[-Awsin(wt  +  y)] 
in  a  Fourier  series  and  retaining  the  first  term  only  (Reference  16). 
The  following  results  are  obtained  for  the  amplitude  A  and  the  phase  y. 


A  = 


-a02  +  y0l  S2  +  a2 (02  ~  0 1 0 3 ) 
0 1 


(58) 


D  =  ^ [A2  ( 1  -  vrr^/mi)2  +  (Ari  +  a/k)2]/(l  +  n2) 


(59) 


Avm2n  -  a (1  -  v  +  n2) 

Y  =  arctan  -  (60) 

Au)2[viri2  -  (1  +  n 2 )  (m  1  +  m2)]  -  nav 


6  =  arctan 


Anm2to2  +  a 


Am2uj2  -  Ak(l  +  n2)  -  n2 


(61) 


where:  0\  =  {mj  +  m2  -  vm2 ) 2  +  n2 (mj  +  m2)2]/(l  +  n2) 
02  =  nkv2  /(I  +  n2) 

03  =  [  (v  -  l)2  +  n2 ] /  (1  +  n2) 


(62) 

(63) 

(64) 


a  =  4yN  R/ir 
v  =  mi'jo2/k  =  ( to /co  1  )  2 


The  solution  (Equation  57)  with  amplitude  A,  expression  (Equation  58), 
exists  when  S/a>/0^  .  From  this  we  can  define  the  coefficient  of  the 
slip  threshold  as: 


6  =  S/a  - 


(65) 
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When  n  =  0,  the  nonlinear  equations  (Equations  58  and  59)  reduce  further 
to  the  very  simple  form  (Reference  17) 


In  these  equations  R  depends  on  the  blade  and  root  geometry.  For  example, 

in  the  earlier  investigations  for  a  twisted  blade  of  length  L,  twist 
1 

angle  a  between  root  and  tip,  dovetail  angle  \ p  and  dovetail  radius  R  , 
it  was  shown  that: 

R  =  Rq  cos  ip/L 

r 

and  that  S  should  be  replaced  by  S  cos  a  . 
b.  Case  2,  Q<Qq,  8  <_  0 

In  some  range  of  the  parameters,  i.e.,  when  8^0,  the  solution 
(Equation  57)  does  not  exist.  This  corresponds  to  the  domain  of  existence 
of  the  linear  solution  for  the  one-mass  system.  In  this  case,  Equations 
51  to  53  become: 

niiXi  +  K  (X 1  +  -  X i )  =  S  cos  ait  (68) 

1  1  1  O) 

X2  *  0  (69) 

The  solution  of  Equation  68  is: 

X  (t)  =  D*  cos  (<*>t  -  8*)  (70) 

with  D*  =  S/k/(l  -v)2  +  n2  (71) 

and  6*  *  arctan  n/ (1  -  v)  (72) 
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c.  Case  3,  Q.<Qq,  B-j  >  0 

Slip  now  occurs  at  the  mass  and  ceases  to  occur  at  the  mass 
m 2>  which  becomes  "locked"  at  high  rotation  speeds.  The  slip  at  the 
mass  m^  can  contribute  significant  amounts  of  damping  if  the  term 
uN^R^  in  Equation  53  can  be  made  to  remain  relatively  small  through 
proper  control  of  the  blade  root  geometry. 


The  equations  now  take  the  form: 


*1  +  *1^  X^  +  *1 

m3X3  +  KljJ<X3 


kJJ  X 3  +  =  s  cos„t 

x !  >  +  (X3  -  X !  )J 


+  K2^  X3  +  X3^  +  uNRj sgnX3  =  0 

The  solution  of  these  equations  is  written  in  the  form: 

X1  =  DjjCOs(ut  +  Yjj  -  <$jj)  and  X^(t )  =  AjjCOS(ojt  +  Yj j ) 


Then,  again  by  the  method  of  harmonic  balance,  it  can  be  shown  that: 


-aj05  +V01+S2  +  ai(05  -  040s) 


f Ajj(1  +  k2/k1)2  +  [Ajjnd  +  k2/k1)  +  o^k^2 
(1  +  n2) 


'  vim3/ml 


Yj-j-  =  arctan 


AXIvlra3h  _  aid  -  vj+  ti2) 

Attu2 [v im3  -  (1  +  n2) (mi  +  m3)  -  naiv; 


Ajj.nni3wz  +  a  i 


6jj  =  arctan 


ATTm3o)2  -  A_tKi  (1  +  n2)  ~  na : 
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vj  -  m1u)2/k1  =  ( 10/10 i  i )  2  (84) 

The  solution  (Equation  75)  exists  when  S/a-,>^  .  From  this  we  define 
the  appropriate  coefficient  g  of  the  slip  threshold  by: 

e  =  S/ a  i  -  /07  (85) 

d.  Case  4,  Q  >  g<0 

In  some  range  of  parameters,  i.e.  B£0,  the  solution  (Equation  75) 
does  not  exist.  It  corresponds  to  the  linear  case,  for  which  the  relevant 
solution  is: 


Xj (t)  =  D*  cos  (wt  -  6*) 

(86) 

with 

Di  =  S/kj/d  -  v  i )  2  +  n2 

(87) 

and 

★ 

6 i  =  arctan  n/(l  -  vj) 

(88) 

Figure  39  illustrates  a  graphical  method  for  constructing  the  solution 
from  these  equations,  for  the  case  n  =  0.  In  this  figure,  note  that: 
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Figure  39  illustrates  clearly  the  regions  of  existence  of  the  nonlinear 
solutions  ( ( A {  f  0)  for  various  values  of  S/ct-j .  When  n  =  0  and  m3  =  0, 
the  solution  further  simplifies  to: 


/l  -  (cti/S)  *  (1  -  Vl)  2 

- -  (90) 

(kj  +  k2)  1  Jc/k !  -  vj  -  V]  (ma/raj)  (k/k2)  (1  -  v  t )  | 


Di 

s 


1  [1  -  v  !  (k/k2 )  (m3/mi  )  ]  2  [  1  -  (aj/SjMl  -  Vj)2]  /alX2 

— \ -  +(— )  (91) 

kj  1  [k/ki  -  vx  -  vj (ma/mj )  (k/k 2 )  (1  -  vj ) ] 2  \S  / 


2  2  ) 

where  it  is  recognized  that  =  w  /w^  and  w-ji  =  /k^Tm^-  For  most 
practical  cases  R,  «  1  in  these  equations  and  for  a  twisted  blade,  we 

*  t 

replace  S  by  S  cos  a.  We  now  have  a  formal  solution  for  the  cases 
n  <  f2Q  and  _>  r2Q »  and  a  numerical  example  will  be  discussed  after  a 
review  of  the  static  blade  behavior.  Finally,  note  that  the  problem 
addressed  here  is  quite  similar  to  that  discussed  by  Williams  and 
Earles  (Reference  3). 

3.  ANALYSIS  OF  QUASI-STATIC  BEHAVIOR 

Referring  to  Figure  40,  the  radial  movement  Y  of  the  sub-platform 

2 

under  the  action  of  the  centrifugal  load  WR^Q  due  to  the  outboard  part 
of  the  blade  will  be  less  than  the  gap  A  up  to  the  speed  ft  where: 

Qo  =  Ada/WrD  (92) 


where  R^  is  the  average  radius  of  the  blade  relative  to  the  rotation  axis 
of  the  disk.  When  Q>i2Q,  the  springs  k^  will  come  into  play  and  provide 
the  normal  load  on  the  mass  m^.  It  is  easily  seen  on  the  basis  of 
static  equilibrium  that: 


WRd(si2  -  nj) 

1  +  kD/2  kp 


(93) 
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which  means  that  kp  >  kp  if  a  significant  reduction  of  the  centrifugal 
loads  on  the  blade  is  to  occur,  as  is  necessary.  As  an  illustration, 
consider  the  blade  for  which  the  static  extensional  stiffness  kp  of  the 
root  below  is  provided  by  a  uniform  segment  of  cross-sectional  area  Ap 
and  length  £p.  Then: 

kD  -  E  V'd  l9,) 

Similarly,  if  each  spring  kp  is  considered  to  be  represented,  as  an 
approximation,  by  a  cantilever  beam  of  length  £  ,  thickness  hp  and 
breadth  bp,  then: 

kp  =  3  E  bph|/12£|  (95) 

For  &<Qq,  the  normal  load  N  on  the  mass  m^  is  given  by: 

N  =  W  (96) 

Table. 2  gives  a  calculator  program  to  predict  the  blade  response 
with  slip  at  the  platform. 


Figure  40.  Quasi-Static  Model 
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TABLE  2 

HP  67  PROGRAM  FOR  BLADE  WITH  SLIP  AT  PLATFORM 


Input ;  Put  a)  in  X  register  (Hz)  ;  press  A,  B,  C,  D  in  sequence. 

Manual  Storage:  Storage  Location  I  -  n  i 

Storage  Location  0  -  mi,  (Kg) 

Storage  Location  1  -  K^,  (N/m) 

Storage  Location  2  -  K2,  (N/m) 

Storage  Location  4  -  c»i/S 
Storage  Location  5  -  S  cos  a 


Output ;  Label  A  - 
Label  B  - 
Label  C  - 
Label  D  - 
Internal  Storage : 


Program; 


|  X]_  |  linear  (a^/S  =  °°)  (Equation  87) 

|  Xj_  |  linear  (a^/S  =  0)  (Equation  71) 
j X^ |  slip  (Equation  90) 

| X2 1  slip  (Equation  91) 

Storage  Location  6  w2  (rad/sec) 

Storage  Location  7  K 

Storage  Location  8  (w/tog)2 

Storage  Location  A  l-fa^/S)2  (l-u2/“i2)  2 


fCL  PRGM 

X 

RCL  6 

X 

1 

RCL  5 

f  LBLA 

gx2 

RCL  0 

hRTN 

+ 

X 

2 

RCL  E 

X 

f LBL  C 

STO  A 

RCL  1 

X 

+ 

CHS 

RCL  1 

RCL  7 

~ 

hir 

f/5 

RCL  7 

RCL  0 

RCL  1 

hRTN 

X 

h  1/x 

+ 

■5 

T 

fLBL  D 

gx2 

RCL  5 

gx2 

h  1/x 

RCL  8 

RCL  7 

STO  6 

X 

STO  E 

RCL  6 

- 

RCL  5 

RCL  0 

hRTN 

RCL  1 

X 

gx2 

X 

X 

f LBL  B 

hRC  I 

STO  8 

h  1/x 

RCL  1 

CHS 

RCL  1 

X 

CHS 

RCL  A 

RCL  1 

h  1/x 

gx2 

1 

X 

RCL  2 

+ 

RCL  2 

RCL  E 

+ 

STO  B 

gx2 

h  1/x 

+ 

RCL  4 

RCL  4 

RCL  B 

STO  E 

+ 

f/* 

X 

gx2 

f  /x 

RCL  1 

h  1/x 

h  1/x 

gx2 

+ 

X 

hRC  I 

STO  7 

RCL  5 

CHS 

f  /x 

hRTN 
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4.  NUMERICAL  ILLUSTRATION 

The  equations  derived  in  this  section  permit  one  to  predict  the 
effects  of  rotational  speed  on  the  dynamic  response  of  any  blade  provided 
that  one  knows  the  modal  masses  and  stiffnesses,  and  the  exciting  forces. 
Unfortunately,  this  information  is  rarely  available  at  the  time  it  is 
needed.  This  is  especially  true  of  the  exciting  forces.  As  an  example, 
the  work  of  Hansen,  Meyer  and  Manson  (Reference  2)  represents  a  very 
interesting  early  investigation  of  blade  damping  as  a  function  of 
rotation  speed.  The  blade  geometry  data  (Reference  2)  indicates  that 
they  used  a  rectangular  section  blade  of  length  50.8mm  (2  inches),  breadth 
15.3mm  (0.602  inches),  width  1.83mm  (0.072  inches)  at  the  tip  and  3.05mm 
(0.120  inches)  at  the  root.  The  root  was  cylindrical,  of  diameter  7.87mm 
(0.310  inches).  In  the  tests,  the  blade  was  attached,  alone,  in  a  33  cm 
(13  inch)  diameter  disk  and  excited  by  impacting  it  with  a  falling  steel 
ball.  Although  the  results  are  very  interesting  and  useful,  it  is  diffi¬ 
cult  to  obtain  accurate  excitation  force  data  from  this  source. 


Ewins  (Reference  20)  recently  described  a  test  system  in  which  air 
jets  were  used  as  a  means  of  exciting  the  blades.  The  test  results  indi¬ 
cate  good  qualitative  agreement  between  linear  analysis  and  experiment  for 
several  tuned  and  untuned  blades  in  a  flexible  disk,  but  no  estimate  of 
the  exciting  forces  is  given  directly.  Other  investigators  (References  14, 
21)  are  equally  uninformative  as  far  as  this  aspect  of  the  problem  is  con¬ 
cerned.  References  22-24  do  seem  to  address  the  question  to  some  extent. 


The  cyclic  forces  acting  on  a  rotating  blade-disk  system  arise  as  the 
blades  cut  through  a  quasi-stationary  airflow  pattern  generated  by  the 
fixed  blades  (vanes)  ahead  of  them.  The  stationary  pressure  field  along  a 
circular  path  through  the  center  of  each  stationary  vane  possesses  a  mini¬ 
mum  between  each  vane  and  a  maximum  at  each  vane  station,  as  illustrated 
in  Figure  41.  If  the  number  of  fixed  vanes  is  n,  then  the  rate  of  repe¬ 
tition  of  the  pressure  pulses  crossing  each  rotating  blade  is  nft/60  Hz,  if 
n  is  the  rotation  speed  in  rpm.  A  Fourier  expansion  of  this  repeating 
pulse  then  gives: 


cos (mnfit/60 ) 


S(t)  =  S  (V ,  ) 


(97) 


nxei> 

VANES 


ROTATING 

BLADES 


I 


FIXED 

VANES 


*  i  \ 

^ — - 
1  I 


along  blades 

PRESSURE  PATTERNS 


P^“  *  - * 

N  AIRSTREAM 


Figure  41.  Flow  Induced  Loads 
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where  =  1  if  the  pulses  are  sharp.  The  amplitude  S(V,ft)  depends  on 
the  mean  velocity  V  of  the  airflow  through  the  stage,  which  in  turn 
depends  on  the  power  setting  of  the  engine,  as  well  as  the  rotation 
speed  ft.  It  is  not  easy  to  determine  S (V,ft)  analytically  and  little 
experimental  data  from  industrial  sources  seem  to  have  been  published. 

So  we  shall  consider  only  a  few  "typical"  cases  in  this  section  in  order 
to  illustrate  the  effect  of  this  important  parameter  on  the  blade  response. 
Each  term  of  Equation  47  gives  rise  to  a  possible  excitation  of  a  blade 
mode,  and  a  typical  Campbell  diagram  is  constructed  by  plotting  the 
frequencies  mnft/60  against  ft  for  various  values  of  m,  along  with  the 
blade  resonant  frequencies  ( f ■  =  1,2)  as  a  function  of  ft,  as  illustrated 
in  Figure  42.  ■  As  is  seen,  the  fundamental  mode  is  excited  by  the  blade 
passage  excitation  when  ft  is  such  that  f^  =  1  nft/60  (nth  engine  order, 
point  A),  by  the  second  harmonic  (2nth  engine  order)  when  f^  =  2  nft/60 
(point  B)  and  so  on.  The  Campbell  diagram  is  a  useful  means  for  esti¬ 
mating  where  the  vibration  problems  are  likely  to  occur,  but  it  gives 
no  clue  as  to  excitation  force  magnitudes. 

We  shall  now  examine  a  specific  blade  geometry,  considered  in 
previous  investigations  (References  17,  18)  and  representing  a  typical 
low  pressure  compressor  blade.  The  main  dimensions  of  this  blade  are: 

L  =  20  cm 

W  =  0.20  kg  (0.44  lb) 

The  other  relevant  magnitudes  are  given  in  the  nomenclature,  as 
specific  numbers  associated  with  the  relevant  symbols.  From  these 
numbers  we  see  from  Equations  94,  95,  92,  93,  and  96  in  succession  that: 

kQ  =  2.10  x  1 09N/m  (1.2  x  107  lb/in) 
kp  =  2.1  x  107N/m  (1.2  x  105  lb/in) 

A  =  0.057  mm  (0.00225  ins) 

N]  =  1.856  x  10"6  (ft2  -  ft2)  lbf 
N  =  8.342  x  10'4ft2Newton(l .875  x  10'4ft2  lbf) 
with  ftQ  =  12000  rpm 
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From  Equations  68,  69,  90,  and  91  we  can  then  calculate  the  response 

of  the  blade  for  the  cases  iKSr  and  .  Some  results  are  shown  in 

o  —  o 

Figures  43  to  46.  From  these  we  can  then  determine  the  apparent  modal 

r 

damping  ,  defined  by  the  "half  power  oandwidtn'  method,  and  plot  against 
a/S  (for  il<S3Q)  or  c^/S  (for  :i>Oo).  The  results  "or  this  particular  blade 
are  shown  in  Figure  47.  It  is  seen  that,  for  a- 5,  ,  n,  increases  rapidly 
as  u/S  falls  and  reaches  a  very  high  maximum  value  before  finally  falling 
to  zero  as  a/S  -*•  0,  i.e.  ->  0.  On  the  other  hand  for  ,  the  damping 
rapidly  increases  as  a^/S  rises  above  a  threshold  level  (below  which  a 
new  peak  occurs  at  lower  frequency  and  with  low  damping),  reaches  a  peak 
of  about  0.20,  and  then  drops  more  slowly  as  ou/$  increases  further. 

From  the  values  of  N  and  given  earlier,  we  can  then  calculate 
a/S  and  a,/S  for  any  assumed  values  of  S.  We  shall  consider  two  cases 

'  r  o  o 


namely  (i )  3  =  22.24  N  (5  Ibf)  and  (ii)  S  =  1.69  x  10 


-6  .  3/2 


(=  3.8  x  10”D  Q,  J/  Ibf),  with  Q  in  rpm.  The  first  represents  a  constant 
force,  as  in  Hanson,  Meyer  and  Manson's  experiments  (Reference  2)  and 
the  second  represents  a  more  likely  situation  in  wnich  S  increases  as  i] 
increases,  e.g.  S  =  22.24  N  (5  Ibf)  at  12000  rpm,  ;;'l.'4  M  (7.0  !bf,.  at 
15000  rpm,  and  5.96  N  (1.34  Ibf)  at  5000  rpm.  Tables  2  and  3  show 

t 

typical  calculated  values  of  a/S,  o^/S,  and  .  The  appropriate 
expressions  for  .,/S  and  ct^/S  are'. 


( i )  S  =  22.2  N(5  Ibf) 


4uN  R 


4  x  0.15  x  1 .875  x  10~4T2  x  1 


S  /rS  Lcosa  cosil/ 


20  x  n  x  5  x  0. 724  x  0.5 


9.89  x  lo'w 


4  l,  N, 


4  x  0.15  x  1.856  x  10“b  -  '/) 


n  x  5  x  0 . i 24 


9.79  x  10 


cos 


RESPONSE  |X,/S  AND  X2/S|  -  m/N  ^  RESPONSE  D/S  OR  U /s|  -  m/N 
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Figure  48  shows  the  variation  of  n-j  with  Q  for  these  cases.  It  is 
seen  that  high  damping  occurs  at  high  speed,  as  expected.  Obviously,  the 
particular  numerical  values  can  be  changed  to  represent  other  conditions, 
but  the  approach  and  general  behavior  remain  the  same.  One  other  possi¬ 
bility  can  be  examined,  namely  a  blade  of  the  type  having  a  hinged  root 
with  mating  surfaces  now  parallel  to  the  plane  of  rotation  of  the  root, 
and  discussed  also  by  Hanson,  Meyer,  and  Manson  (Reference  2).  For  this 
type  of  blade,  N-j  remains  just  about  constant  at  all  speeds,  white  the 
force  S  changes,  usually  increasing  as  Q  increases.  This  would  account 
for  the  high  damping  achieved  by  this  type  of  blade,  since  a,  would  be 

t  ' 

nearly  constant  ana  a,/S  would  change  in  such  a  way  that  n-j  would  increase 
as  9,  increases. 


S  =  22  24  Newtons 

l  O-O-  S  ;I69  Xl0’’n5'2  Newton* 


ROTATIONAL  SPEEO  it  KRPM 


Figure  48.  Modal  Damping  versus  Rotational  Speed 
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SECTION  IV 

HARMONIC  RESPONSE  OF  A  DAMPED  TWO-DEGREE  OF  FREEDOM 
SYSTEM  WITH  GAPS 


1.  INTRODUCTORY  REMARKS 

In  this  section,  a  two  mass  mechanical  system  excited  by  an  external 
harmonic  force  is  considered.  One  of  the  masses  is  assumed  to  move  in 
the  space  between  two  elastic  non-inertial  buffers.  The  gaps  introduce 
a  high  degree  of  nonlinearity  into  the  system.  The  motion  of  the  two 
masses  is  damped  by  dry  friction  which  also  brings  nonlinear  effects  and 
hysteretic  damping  of  the  elastic  members.  An  approximate  steady  state 
solution  is  found  for  the  harmonic  response  case.  A  graphical  method 
is  adapted  for  construction  of  response  (amplitude  versus  frequency) 
curves.  Numerical  examples  show  the  influence  of  various  parameters, 
such  as  amplitude  of  the  exciting  force,  dry  friction,  nysteretic  damping, 
fluid  pumping,  and  elasticity  of  the  buffers  on  the  amplitude  of  the 
response.  This  section  summarizes  Reference  16. 

It  is  generally  very  difficult  to  predict  the  dynamic  oenavior  of  a 
mechanical  assembly  having  gaps  between  the  members,  whether  introduced 
for  specific  reasons  or  resulting  from  wear.  In  fact,  the  motion  of 
systems  with  gaps  is  highly  nonlinear  and  is  very  sensitive  to  the  gap 
dimensions  and  the  values  of  the  excitations  (Reference  35). 

The  present  section  considers  the  response  of  a  damped  two  degree 
of  freedom  system  with  gaps,  excited  by  a  harmonically  varying  force. 

The  model  of  the  system  is  based  on  a  modelization  of  blades  in  a  disc. 
Damping  in  the  system  is  introduced  by  dry  friction  between  one  of  the 
masses  and  the  supporting  environment,  hysteretic  damping  accompanying 
motion  of  the  elastic  members  and  fluid  pumping  in  gaps. 

A  combined  analytical-graphical  method  is  adopted  for  determining 
the  response  curves,  the  amplitude-frequency  relationship.  This  allows 
one  to  determine  the  steady  state  response  of  the  system  for  various 
values  of  the  important  parameters,  particularly  gaps,  damping,  friction, 
and  excitation. 
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2.  MODEL  OP  THE  SYSTEM 

The  idealized  model  of  the  system  is  illustrated  in  Figure  49.  The 
system  consists  of  two  rigid  bodies  of  mass  M,  and  M.,  ioined  by  a  plane 
linear  spring  of  rigidity  and  two  linear  buffer  springs  of  rigidity  k^. 
Assumptions  made  in  the  analysis  are:  (i)  the  system  is  symmetric  and 
plane,  (ii)  the  masses  of  the  buffers  are  negligible,  and  (iii)  the  body 
M^  moves  horizontally  only,  while  the  body  M,  performs  plane  motion. 

The  external  exciting  force  S  cos  (wt  -  y)  is  applied  to  the  body  . 
During  the  motion  of  the  body  M^  a  dry  friction  force  arises,  modelled 
by  Coulomb's  law  (F  =  pN,  where  u  is  the  coefficient  of  friction  and  N 
the  normal  force).  When  the  amplitude  of  motion  of  the  body  M^  is  greater 
than  the  gap  A,  the  springs  k^  exert  a  force.  The  remaining  damping 
terms  in  the  system  are  linear.  For  small  motions  of  the  bodies,  the 
equations  of  motion  of  the  system  are: 

M,  (R<Jj+X)  +  ($+  5)  ■*  M=S  COS  (wt-Y)  (98) 

I  CO  K  S 

M2X  +  “x*  +  [  ~  X  +  k2(x-A)  ]  0 

-k^  +  pN  sgn  (X)  =  0,  (99) 

where 

'P  =  |  [1+sgn  ( | X 1  -  A)] 


ax  represents  the  viscous  damping  coefficient  (e.g.,  fluid  resistance), 
a^,  01^  represent  the  hysteretic  damping  coefficients  (e.g.,  internal 
friction  in  the  springs).  We  will  confine  attention  to  the  effects  of  gaps 
and  damping  on  the  response  of  the  system. 


60 


AFWAl-TR-80-4003 


Figure  49.  Model  of  the  System 

3.  STEADY-STATE  SOLUTION  WITH  SLIP 

We  shall  analyze  the  excited  motion  only,  i.e.: 

X  =  Acosuit, 

where  A  represents  the  amplitude  of  the  excited  vibration.  I 
the  frequency  of  the  exciting  force. 

Putting  Equation  100  into  Equation  98  we  obtain 

a  . 

M1R^  +  ^  +  =  Scos  (wt-y)  + 

+  M^Aw^  coswt  +  ^  sin  cut 

The  steady-state  solution  for  Equation  101  is 

'i'(t)  =  p  ^  (l-  cos  (wt-y)  +  sin(ut-Y 

kp*  1  1 

,  A 1  T  /  ,  M1R  2 \  .  ajjjh  .  .I  A 

+  \  1-  — r~  w  1  coSwt  +  --t—  smut  -  —  rnsi 


(IPO) 

t  depends  on 


(101 

»] 


(in? 
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where 


Putting  Equations  100  and  102  into  Equation  99  we  get: 


b  1,  ‘~ 

r  ?  Ki  *1 

A  cosut  1  k~ <l>*  -  M ~u  +n - “y 

(ti. 

VR 

4>* 

-  -  A  si  nut  £  ctxw+  ctxj1  y*  + 

kla^h  1 

R4>*  * 

k^Ay*  -  uN  sgn  (A  sinut)j  - 

k1<P* 

[<kr 

M1RuJ2)cos  lut-y)  +  Sin  (ut-y)  j 

=  0 

(10 

where: 

T*  =  [l-sgn(  |Acosutj-A)]/2  '10 

The  values  of  the  function  y*  are  given  in  Table  4. 

We  develop  the  function  y*,  y*  coswt,  y*  sinwt  and  sgn  (A  sinut)  in 
Fourier  series.  Limiting  attention  to  the  first  approximation  we  have: 

2t r  2ji 

i'*  =  '  [  coswt  y*cosut  dt  +  si  nut  J u  y+  si  nut  dtj  =0  (10 

0  0 

o 

Note  that  when  ^,.=0,  >*=1  -M^Ruf/k^ . 
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''co Sait  =  ^  j^cosuit  J  w  f*  cos2uit  dt  +  sinuit  4**  sinuitcosuit  dtj 


=  ^  (2uit-j  +  sin  2uit^)cosuit  =  C(A)cosuit 


4'*  sin  wt 


^cosuitj"  w  41*  sinuit  cosuit  dt  +  sinuit  J u  4'*  sin2uit  dtj 


( 2oit i  -  sin  2oot i )  sinuit  e  C*  (A)  sinuit 


sgn  (Asinuit)  =  sgn  A  ^  ^  cosuit  J  w  sgn  (sinuit)cosuitdt 


sinuit  J  w  sgn  (sinuit)  sinuit  dtj 


sgn  A  sinuit 


Putting  Equation  105  to  Equation  108  into  Equation  103  and  comparing  the 
coefficients  of  the  trigonometric  functions  we  obtain  the  basic  equation 
for  the  resonance  curve  A  =  A(ui).  (The  same  result  is  obtained  if  one 
assumes  the  solution  in  the  form  x  =  A  cosuit,  4*  =  Bcosuit  +  D  sinuit,  where 
A,  B,  D  are  functions  of  ui). 

A%  (ui)  +  [A<£2  (ui)  +  BStJi*  (ui)  sgn  A]2  =  S2  (109) 
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and  the  expression  for  the  phase  angle: 


where 


(|s<))*sgnA  +  A<j>2J(k, -M,Rw  )  +  Aa  h«^j" 

Y  =  arctan  t-t - , - o'—* - 

^S4>*sgnA  +  A4>2^ ot  h -A ( k ]  -M1  Roj 


(110) 


(111) 


(112) 


(113) 


The  difficulty  of  the  problem  arises  because  the  expression 
(Equation  109)  is  nonlinear  with  regard  to  A,  due  to  functions  (sgn  A) 
and  t^  =  1 1  (A).  The  last  function,  which  is  included  in  C  (A)  and  C*  (A), 
depends  also  on  the  gap  A  (Equation  105).  The  functions  sin  2wt^+2ojt^ 
and  2u>t|-sin  2u>t^  are  presented  in  Figure  50.  As  we  see  in  Figure  50 
the  function  (sin  2u)t^+2wt^)  is  practically  constant  for  A>2A  and  equals 
about  3.1. 

We  shall  seek  for  the  explicit  solution  A=A(w)  of  Equation  106 
supposing  that  the  functions  C(A)  and  C*(A)  (Equations  103  and  109)  are 
constant.  Then  we  shall  introduce  the  functional  dependences  C=C(A), 

C*=C*  (A),  applying  a  graphical  method. 

If  the  C  and  C*  are  supposedly  constant  we  obtain  from  Equation  62: 

ft  <J>  *  <t>  ?  -  J $  i +  f  ~  /  ? 2  $  ^ 

A(u)  =  +  S  - - -  (114) 

MU  !  +4=  2  ) 


which  is  valid  for  -  1 . 
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Figure  50.  Graph  of  t-j  =  l/w  Arc  Cos  A/A 


As  we  see,  for  large  values  of  the  excitation  amplitude  S,  there 
exists  one  positive  value  for  A  and  one  negative  of  the  same  value.  For 
small  values  of  excitation,  motion  of  this  kind  is  not  possible,  and 
only  ip  ( t )  exists  -  linear  motion  of  the  mass  . 

The  expression  (Equation  114)  represents  the  relation  A=A(w),  i.e., 
the  relationship  between  the  amplitude  of  the  excited  vibration  and  the 
frequency  of  the  exciting  force  in  the  simplified  case  C=const, 

C*=const.  A  graphical  method  of  solving  Equation  109  with  C  and  C*  not 
constant  will  be  explained  for  the  particular  case  of  a  system  with  dry 
friction  only  in  Section  IV. 5. 


4.  STEADY-STATE  SOLUTION  WITHOUT  SLIP 

Finally,  for  the  case  of  small  exciting  force  amplitude  S,  that  is 
S<4uN(f*/ir,  no  slip  of  the  mass  occurs,  unless  >  pN.  Then  the 
motion  is  linear  and  involves  ^  alone,  so  that  the  equation  of  motion 
of  the  system  becomes 


M  i  R<|j 


'fh 


iii  +  ki^=  Scos(wt-y) 


(115) 
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This  can  be  solved  directly  for  i>4t).  The  graph  of  (where  <i>*  is  an 
amplitude  of  the  response  of  \p)  versus  w  is  then  a  classical  response 
curve  up  to  the  point  where  <J;*>  pN/k^,  at  which  point  slip  occurs,  but  a 
sustained  solution  with  slip  does  not  exist  (because  S<4uN4>*/7r) .  There¬ 
fore,  the  response  under  these  conditions  is  "cut-off“  at  the  level 
•y  =  pN/k-| ,  i.e.  amplitudes  are  limited.  Only  for  S>4pN^*/tt  can 
infinite  amplitudes  be  sustained  at  resonance,  as  in  Equation  114,  for 
zero  values  of  the  other  damping  coefficients. 


5.  SOME  SPECIAL  CASES 


Ihe_«se 
without  damping. 


=q^=a^=0  -  Dry  friction  only. 
The  expression  (Equation  114) 


Let  us  consider  the  case 
has  then  the  form: 


A(u) 


(116) 


where 


1  - 


-  1- 


.Ml+M2  kzC  kzC 

*~W2( - n -  +  ~T7-)  + 


MjR 


M 

.,2 


M; 


M 


(we  omit  the  sign  minus  in  front  of  the  expression  (Equation  116);  of 
course  there  are  two  solutions:  positive  and  negative). 

The  function  A(uj)  exists  only  for 


0 


M,  R 


f  1+  77  I  if 


M  1 

<7*  i1' 


M,R 


1+  l) 


if  ay l 


2  *  i  ] 

and  if  (1  f  -;)  then  it  takes  the  value  zero. 

M^R 
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The  function  A(w),  (Equation  116),  tends  to  infinity  when  oj  -*  or 

tu  -*■  n>2  >  where 


For  C=C  ,  where  (Figure  50) 
m3  x 

J  (2u.t,tsin2u,t,)max  -  3J  -  0.987  (1  18) 


we  have 


For  u>  =  0,  A(0)=S/fTp2Vk,,C  which  -+  if  C  ->  0. 

Finally  the  function  A(w)  (Equation  116)  to*"  some  cases  has  shapes 
presented  in  Figures  51,  52,  and  53. 

As  one  sees,  the  presence  of  the  dry  friction  does  not  eliminate  the 
resonance  peaks.  The  curves  presented  in  Figures  51,  52,  and  53  are 
plotted  for  the  values  C  =  const.  In  reality  the  value  of  C  varies: 

0  -  C  ■'  Cmax  =  0.987,  when  A  varies  in  the  range  A  ^  A  ■_  2A.  Then  for 
A  A  the  curve  A(u)  with  C  =  0  is  applicable.  For  A  2a,  the  curve  with 


of  Response  Curve  for  the  Nonlinear 
Case 
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C  =  Cmax  is  used.  For  A  <  A  <  2A  one  should  plot  a  family  of  curves  A(u) 
with  a  parameter  C  varying  0  <  C  <  Cmax  and  find  the  corresponding  points 
at  intersections  of  those  curves  and  lines  A  =  R.A,  1  <  R.  <  2, 

l  i 

i  =  1,2,3,...,  the  values  R.  being  taken  from  Figure  50. 

Figure  54  presents  the  method  of  plotting  the  real  curve  A  (w) 
(Equation  116).  The  shape  of  the  A(ta)  curve  depends  on  the  value  of  A. 

As  can  be  seen  in  Figure  54,  the  resonance  peaks  of  the  A(w)  curve 
(the  case  without  damping)  always  exist  (two  or  one  only).  The  dry 
friction  does  not  cancel  the  resonance,  it  reduces  only  the  amplitude  A. 

The  opposite  effect  is  brought  into  play  by  the  amplitude  of  the  exciting 
force  S.  The  existence  of  the  gap  A  changes  qualitatively  the  A(co)  curve. 
It  brings  a  highly  nonlinear  effect.  For  limiting  cases  A  =  0  and  A  =  °°, 
the  response  curves  correspond  to  the  cases  C  =  0,  C  =  1  (Figures  52  -  54) 
and  then  there  is  no  nonlinear  effect.  In  other  cases,  the  nonlinear 
effect  on  A(w)  curves  is  more  distinct  if  A  is  greater. 

Numerical  Example  No.  1 

Let  us  consider  the  example  with  the  following  data 

M-,  =  1 .827  kg,  =  15.3  kg, 

k1  =  14272.7  N/rad,  R  =  9.09  10~3m. 

Four  values  of  B:  6  =  0,  0.2,  1,  1.1  and  three  values  of  k^  were  considered 

k2  =  1 06N/m;  107N/m;  108N/m.  (120) 

The  calculations  were  made  on  an  SR  52  calculator. 

The  resonance  curves  (Equation  116)  have  been  plotted  in  coordinates 

(ro,AM/S).  In  Figures  55-57  we  have  the  family  of  curves  (Equation  116) 

with  different  constant  parameters  (C/0  <  C  <  C  )  for  the  case  n  =  0. 

They  correspond  to  the  values  A/A  =  1 ;  1 . 1  ;  1 . 2 ; /2;2 .  One  can  see  the 

normal  shape  of  the  linear  resonance  curves  for  a  two-degree  of  freedom 

6  7 

system.  For  the  cases  k-  =  10  N/m  and  k„  =  10  N/m  the  first  mode  corre- 

C.  L  O 

sponds  to  the  motion  of  the  mass  M^.  For  the  case  k^  =  10  N/m  the  first 
mode  is  that  of  the  mass  .  The  second  resonance  zone  in  the  cases 
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Figure  56.  Response  Curves  for  C  =  co 
(k2  =  107  N/m,  3=0) 
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k^  =  lO^N/m  and  k^  =  10^N/m  for  different  C  is  nearly  the  same,  so  in  this 
region  the  gap  nonlinear  effect  will  be  small.  Such  low  values  of  k£  do 
not  alter  much  the  values  and  (Equation  117  and  Figures  51  and  52). 

Then  the  gap  A  is  introduced  and  the  curves  corresponding  to  non¬ 
linear  cases  are  plotted,  applying  the  method  explained  previously.  For 
every  value  of  k^  two  cases  of  the  values  of  the  gap  A  were  treated 
(Figures  58-60).  To  show  the  influence  of  the  exciting  amplitude  the 
previous  curves  (Figures  58a  -  60a)  are  plotted  on  the  plane  (oj.A) 

(Figures  58b  -  60b),  where  the  scale  of  A  is  expressed  by  the  value  of 
the  gap  A.  As  can  be  seen,  a  smaller  value  of  S  causes  a  reduction  of 
the  amplitude  A  and  narrows  the  resonance  zone.  It  cannot,  however, 
cancel  completely  the  resonance  peak.  For  the  smaller  gaps  (Equation  1) 
the  nonlinear  effect,  which  obviously  has  a  "hard"  character,  is  less 
exhibited.  The  gap's  nonlinear  effect  lies  in  the  appearance  of  non¬ 
uniqueness  of  the  amplitude  A  for  certain  regions  of  w.  This  effect  is 
especially  evident  for  large  values  of  k^  (Figure  60).  One  should  expect 
that  the  middle  value  of  the  amplitude  (dashed  curve.  Figures  58b  -  60b) 
gives  an  unstable  solution  (Equation  100). 

The  comparison  between  Figures  58,  59,  and  60  show  the  effect  of  the 
buffer's  rigidity. 

The  curves  for  the  case  with  small  dry  friction  factor  (3  =  0.2) 
are  similar  to  the  undamped  case  (compare  Figures  61  and  60a).  The  effect 
of  the  dry  friction  damping  lies  in  limiting  the  curves  on  the  frequency 
axis.  For  low  and  for  high  frequency  excitation  there  is  no  motion 
(Equation  100)  for  the  mass  (Figures  62-65).  For  the  higher  dry 
friction  some  resonance  zones  are  more  narrow  but  the  nonlinear  effect  of 
gaps  are  more  exhibited  (compare  Figures  61  and  58a).  If  the  dry  friction 
is  sufficiently  high  the  first  resonance  zone  can  be  cancelled  and  the 
effect  of  gap's  nonlinearity  disappears  (Figure  62).  This  effect  does 

O 

not  exist  for  the  case  k^  =  10  N/m  (Figures  63-65)  It  is  obvious,  because 
the  first  resonance  in  this  case  corresponds  to  the  motion  of  the  mass  M-j 
and  there  is  no  possibility  to  cancel  it  by  the  dry  friction  applied  to 
the  mass  M. 
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The  Case  uN 


a  h  =  a.th  =  0.  Fluid  pumping  only. 


Let  us  consider  the  case  without  dry  friction  and  without  hysteretic 
damping.  Only  the  dissipation  by  fluid  pumping  in  gaps  is  present.  The 
response  curve  is  given  by  Equation  116  (for  C  =  const)  with  q  =  0  (we 
omit  the  sign  minus,  of  course,  there  exists  two  solutions:  a  positive 
and  a  negative) : 


For  small  values  of  the  damping  coefficient  a. x  the  response  curve 
(Equation  121)  has  two  maximum  peaks  (for  w  near  and  oj^,  Equation  70), 
Figure  66.  The  position  of  the  peaks  depends  on  the  values  of  cx  and  C 
as  well.  So  the  nonlinear  effect  of  the  gap  can  be  reflected  in  the  A(w) 
curve  in  the  same  manner  as  previously,  plotting  a  family  of  curves  with 
C  as  a  parameter  and  finding  the  intersection  points  with  corresponding 

level s  A  =  R^A,  where  1  <  R^  <  2 ,  i  =  1 ,2 -  For  A  £  A  the  curve 

(Equation  121)  with  C  =  0  is  valid.  For  A  >  2A  the  curve  (Equation  121) 
with  C  =  Cmax  (Equation  118)  is  valid. 


Figure  66.  Resnonse  Curves  for  a  Damped  Case 
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Numerical  Example  No.  2 

Let  us  consider  the  same  example  as  previously,  with  n  =  0,  adding 
the  fluid  damping  coefficient  (four  different  values): 

ax  =  103Ns/m;  1500  Ns/m;  4000  Ns/m;  7000  Ns/m.  (123) 

The  results  of  the  calculations  of  the  response  curves  (Equation  121) 

are  presented  in  Figures  67-72.  The  curves  obtained  applying  the  same 

graphical  method  are  plotted  in  the  coordinates  (to, A),  where  A  is 

expressed  in  terms  of  the  gap's  magnitude  A.  One  can  see  the  influence 

of  the  gap  coefficient  MA/S  on  the  nonlinear  behavior  of  the  A(<o)  curves. 

The  curves  show  the  influence  of  the  amplitude  of  the  exciting  force  S 

on  the  A(io)  relationship.  The  growing  S  causes  an  increase  of  A  and 

some  qualitative  changes  (appearance  of  new  branches  of  A(iu)).  The 

influence  of  the  damping  coefficient  ax  can  be  investigated  by  comparing 

the  curves  corresponding  the  same  value  of  (Figures  68  -  72).  The 

curves  A(<o),  for  small  values  of  ax  do  not  differ  much  from  the  curves  for 

the  previous  case  6=0  (compare  for  instance  Figures  59b  and  68).  The 

greater  values  of  «x  cause  the  disappearance  of  some  A  branches  (compare 

Figures  59  and  69  or  60  and  72),  a  reduction  of  the  amplitude  A  and  a 

disappearance  of  the  resonance  (compare  Figures  58  and  67,  59  and  69, 

60  and  71)  and  a  "softening"  of  the  response  curve  (compare  Figures  69 

and  70).  The  growing  coefficient  a  changes  a  little  the  frequency 

oj-j  ,  of  the  greatest  amplitude  A  (to  smaller  values  of  <o  for  =  10  N/m, 

brings  closer  both  of  them  for  k„  =  10^N/m).  The  increasing  coefficient 

£  8 

ax  cancels  more  efficiently  the  first  resonance  (for  k2  =  10  N/m).  The 
same  value  of  «x  has  more  influence  on  the  shape  of  the  resonance  curve 
if  k2  is  smaller  (compare  Figures  58,  59,  67  and  69). 
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The  Case  N  =  g„  =  ot^  =  0.  Hysteretic  damping  in  the  buffer 
springs  only 

Let  us  consider  the  case  without  the  dry  friction  and  without  the 
fluid  pumping.  Only  dissipation  by  hysteretic  damping  in  the  buffers 
is  present.  The  response  curve  is  given  by  Equation  114  (for  C  =  const, 
C*  =  const)  with  6=0  (we  omit  the  sign  minus): 


where  in  this  case 


A  (u ) 


S 

M 


_ 1 _ 

^<t>  l  (w)  +  <J>2  (w) 


is  the  same  as  previous  (Equation  122)  and 


(124) 


<f>2  (°J) 


a  .  C* 
xh 

M 


M  i  R 


For  small  values  of  damping  coefficient  ax^  the  response  curve  (Equation 
124)  has  two  maximum  peaks  (for  u>  being  near  u>-|  and  see  Equation  117 
and  Figure  66) . 


Numerical  Example  No.  3 

Let  us  consider  example  2  (with  6=0),  adding  the  hysteretic 
damping  coefficient  (three  values): 

axh  =  107N/m;  108N/m;  2  1010N/m. 

The  results  of  the  calculations  of  the  response  curves  (Equation  124)  with 
different  values  of  «x^  are  presented  in  Figures  73-75.  The  curves 
are  plotted  in  the  same  coordinates  (w,  AM/S)  as  previously.  There  can 
be  seen  the  influence  of  the  magnitude  of  the  gap  coefficient  AM/S  on 
the  nonlinear  behavior  of  the  A(oj)  curves.  The  damping  brought  by  the 
buffers  has  nearly  the  same  influence  on  A(w)  relationship  as  the  fluid 
pumping.  However,  the  effect  of  damping  contributed  by  the  coefficient 
is  much  smaller  than  that  of  «x  (compare  the  magnitudes  of  a  ^  and 
<%  which  give  nearly  the  same  effect  of  damping  -  Figures  69  and  73, 

70  and  74,  71  and  75).  For  sufficiently  high  damping,  the  amplitude 
A  of  vibration  does  not  exceed  the  value  of  the  gap  A.  The  nonlinear 
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effect  due  to  the  gap  is  present  in  all  the  cases  (regions  of  frequency 
where  three  values  of  amplitude  are  present,  with  sudden  non-smooth 
changes  of  amplitude  values). 

The  Case  f_a  ^  =  0.  Dry  friction  and  fluid  pumping  only 

Let  us  consider  the  case  without  hysteretic  damping.  The  dissipation 
is  caused  by  dry  friction  and  fluid  pumping  together. 

The  resonance  curve  is  given  by  Equation  114  (for  C  =  const  and 
C*  const) : 


A  (a) ) 


~ /3 2 *  +  i  (l-/?2 <t>l )  +  ^2 


M  ( ( 


) 


(125) 


where  0^  is  given  by  Equation  122  and 

M !  R 

=  1  kj  w2' 

Numerical  Example  No.  4 

We  shall  consider  example  2  with  the  fluid  pumping  coefficient 
(two  values) : 

u  =  2000  Ns/m;  4000  Ns/m. 

The  results  of  calculations  can  be  seen  in  Figures  76  and  77.  The 
common  effect  of  the  dry  friction  and  the  fluid  damping  reduces  the 
resonance  amplitude  and  limits  the  zone  of  frequency  (compare  Figures  61, 
62  and  70,  77).  The  gap's  effect  has  nearly  disappeared. 


4>  = 
2 


“x^* 
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6.  TIME  DEPENDENCE  OF  THE  DEFLECTION  ip 

After  having  found  the  value  of  the  amplitude  A  of  the  solution  for 
given  parameters,  and  for  every  value  of  the  frequency  01,  it  is  possible 
to  obtain  the  time  dependence  of  the  deflection  d>( t ) .  Putting  Equation 
100  into  Equation  102  we  obtain: 

V'  ( t )  =  (k2>?  -  M2to2)  cos  cot  -  ^  f“x  + 

01  h  * 

V  )  sin  ait  -  uN  sign  (sin  ait)  .  (126) 

Taking  into  account  the  values  of  ip*  (Table  4)  for  every  interval  of  time, 
we  can  plot  the  relation  (Equation  126)  in  coordinates  (t,  ip  (t)).  Figure 
78  shows  the  graphical  construction  of  the  Equation  126  for  small  and 
great  values  of  frequency.  As  it  can  be  seen,  due  to  the  nonlinear  effects, 
the  time  dependence  of  the  deflection  \ji(t)  represents  a  complicated  curve. 
Practically,  in  real  systems,  it  will  be  much  more  smooth  -  the  modelization 
accepted  here  is  quite  idealized. 


Figure  78.  Time  Dependence  of  the 

Deflection  of  the  Mass  M1 
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7.  SUMMARY 

The  response  of  a  mechanical  system  to  a  periodic  excitation  is 
strongly  dependent  upon  the  tolerance  between  parts.  In  this  section 
it  has  been  shown  that  the  response  (amplitude  versus  frequency)  is 
very  sensitive  to  the  value  of  the  gaps  -  the  response  curves  change 
not  only  quantitatively  but  qualitatively  as  well  with  variation  of  the 
gap's  magnitude.  The  presence  of  gaps  introduces  nonlinear  effects. 

A  change  of  the  gap's  magnitude  may  increase  the  amplitude  of  the 
response,  which  can  cause  premature  damage  or  failure  of  the  system. 

Another  nonlinear  effect  is  caused  by  dry  friction  in  the  system. 

The  dry  friction  reduces  the  response  amplitudes,  it  cannot  completely 
cancel  the  resonances.  Both  nonlinear  effects  cause  the  nonuniqueness 
of  the  amplitude  of  the  response  for  certain  regions  of  frequency. 

The  external  damping  brings  usual  effects  of  amplitude  reduction. 

A  graphical  method  was  adapted  for  construction  of  the  response  curves 
for  nonlinear  systems.  This  method  is  valid  for  any  nonlinearity  in  the 
system,  but  its  advantages  are  especially  clear  for  discontinuous 
functions. 
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SECTION  V 

EXPERIMENTAL  INVESTIGATIONS 


1.  BLADE  DESCRIPTION 

The  analytical  part  of  this  report  describes  the  fundamental  mode 
dynamic  behavior  of  a  typical  jet  engine  compressor  or  turbine  blade 
having  a  simple  dovetail  root,  with  allowance  for  slip  at  the  root. 

The  purpose  of  the  test  program  was  to  verify,  as  far  as  possible,  the 
results  of  the  analysis.  The  blade  selected  for  the  initial  tests  was  a 
relatively  large  twisted  steel  blade  with  a  simple  root  geometry  as  shown 
in  Figures  79  and  80.  The  blade  was  held  in  a  heavy  fixture  having 
mating  surfaces  to  match  the  contours  of  the  blade  root.  The  centrifugal 
load  was  represented  by  means  of  simple  spring  loaded  wires  applied  at 
the  root  in  such  a  way  as  to  minimize  interference  with  the  blade  root 
motion  during  slip.  This  arrangement,  shown  in  Figure  80,  is  not  ideal 
in  so  far  as  it  does  not  allow  for  the  effect  of  the  centrifugal  load  in 
untwisting  the  blade,  but  it  does  represent  the  root  conditions  far  more 
adequately  than,  for  instance,  clamping  the  root  by  means  of  a  bolt,  which 
would  prevent  slip  and  hence  inhibit  the  very  phenomenon  being  studied. 

The  blade  was  about  20  cm  long. 

2.  TEST  SYSTEM 

For  determining  the  mode  shapes  of  the  blade,  a  digital  test  system 
illustrated  in  Figure  81  was  used.  The  blade  was  solidly  clamped  at  the 
dovetail  and  excited  by  a  small  shaker  through  a  force  gauge  at  the  block. 
The  acceleration  response  was  measured  at  several  points  on  the  blade 
surface  using  a  miniature  accelerometer  and  transfer  functions  measured 
at  several  points  on  the  blade. 

For  the  harmonic  vibration  tests  on  the  blade  with  slip,  the  analog 
test  system  illustrated  in  Figure  82  was  used.  The  Bruel  and  Kjaer  Model 
1014  Beat  Frequency  Oscillator  can  generate  a  harmonically  oscillating 
voltage  of  magnitude  0  to  120  volts  at  any  selected  frequency  from  20  Hz 
to  20  KHz.  The  output  impedance,  measured  in  ohms,  can  be  varied  according 
to  the  type  of  excitation  device  used  in  order  to  match  impedances  as 
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closely  as  possible.  As  the  frequency  is  varied,  the  current  through 
the  output  terminals  to  the  exciter  can  be  controlled  by  means  of  the 
"compressor".  This  is  simply  a  feedback  loop  which  measures  the  voltage 
across  a  fixed  resistor  in  the  output  circuit,  and  uses  it  to  control 
the  output  current.  A  meter  gives  the  output  voltage  across  the  terminals. 

The  output  current  of  the  oscillator  was  fed  directly  to  a  magnetic 
transducer  (Electro  Model  3030-HTB).  The  impedance  setting  of  the 
oscillator  was  set  at  600  ohms  in  order  to  best  match  the  impedance  of 
the  transducer,  thereby  minimizing  distortion  of  the  output  signal.  The 
transducer  consists  essentially  of  a  magnetized  iron  rod  with  many  coils 
of  fine  insulated  wire  surrounding  it,  through  which  the  oscillating 
current  flows.  This  current  produces  an  oscillating  magnetic  field 
which  modulates  the  steady  magnetic  field  of  the  rod,  and  hence  produces 
an  oscillating  harmonic  force  on  any  iron  object  placed  nearby.  The 
magnitude  of  this  force  depends  on  the  amplitude  of  the  oscillating  on 
current  and  on  the  gap  between  the  end  of  the  magnetized  rod  and  the 
object  being  excited.  A  typical  transducer  section  is  shown  in 
Figure  83. 

The  waveform  of  the  input  signal  to  the  transducer  was  monitored 
on  a  dual  -  beam  oscilloscope  (Ballantine  Model  1066S).  The  frequency 
of  the  signal  was  measured  by  a  digital  frequency  meter  (Hewlett  Packard 
Model  5216A).  The  pickup  system  used  a  miniature  accelerometer  (Endevco 
Type  22)  weighing  about  0.2  grams.  The  acceleration  amplitude  of  the 
blade  at  the  point  where  the  accelerometer  was  attached,  with  cyano¬ 
acrylate  adhesive,  was  detected  by  the  quartz  crystal  in  the  accelero¬ 
meter,  unlike  the  transducer,  is  an  extremely  high  impedance  device  so 
that  a  special  high  input  impedance  amplifier  (MB  N-400  Zero  Drive 
amplifier)  had  to  be  used  to  amplify  the  signal  before  it  could  be  read 
off  the  Vacuum  Tube  Voltmeter  (Hewlett  Packard  Model  3400A)  and  monitored 
on  the  other  channel  of  the  dual  beam  oscilloscope.  A  "line-driver" 

(MB  Type  9402215)  was  used  to  reduce  loss  of  signal  in  a  relatively 
long  cable.  Because  of  the  high  system  impedance,  current  levels  in  the 
cables  were  very  low  and  care  had  to  be  taken  to  avoid  excessive  loss  of 
signal  or  pickup  of  stray  signals. 
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Figure  83.  Transducer  Section  (Sketch) 

3.  MEASUREMENT  OF  MODE  SHAPES 

For  the  first  two  modes  of  the  compressor  blade,  the  transfer 
functions  |  X/S|  were  measured,  as  a  function  of  frequency,  using  a  small 
shaker  acting  through  an  impedance  head  for  excitation  and  an  accelero¬ 
meter  for  pickup  as  in  Figure  81.  Results  were  plotted  for  points  (0,0), 
(0,1)  ---  as  illustrated  in  Figure  80.  Figure  84  shows  a  typical  response 
Figures  85  to  89  show  more  detailed  plots  near  each  resonance.  The  system 
was  not  calibrated.  The  product  (X/S)on  was  then  calculated  for  each 
mode,  as  shown  in  these  figures.  Note  that  ( X/S )  is  drawn  on  a  decibel 
scale,  so  that  one  must  convert  back  to  a  linear  scale  in  performing  the 
calculation,  as  follows: 


|  linear  log,'1  (  ) 


For  example,  for  point  (0,0),  mode  1 ,  <r  =  39db  and  =  0.00287. 


'°s;J  (I) 


x  0.00287  =  0.256 


The  summarized  modal  data  is  given  in  Table  5.  These  modes  are  plotted  in 
Figures  90  and  91.  The  modes  can  now  be  used  to  calculate  the  equivalent 
masses  and  stiffnesses  of  the  discrete  model  of  the  blade. 
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Figure  87.  Resonance  Peaks  for  Point  0,2  Figure  88.  Resonance  Peaks  for  Point  3,1 
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Figure  89,  Resonance  Peaks  for  Point  3,2 


TABLE  5 

MODE  SHAPE  DATA 


Point 

(X/S) n i 

(X/S) n  2 

l 

<f>  2 

0,0 

+0,256 

+0,157 

+1,000 

+1,000 

0,1 

+0,234 

+0,117 

+0,914 

+0,745 

0,2 

+0,214 

+0,086 

+0,836 

+0,548 

1,0 

+0,189 

+0,039 

+0,738 

+  0,248 

1,0 

+0,148 

+0,004 

+0,578 

+0,025 

1,2 

+0,147 

-0,023 

+0,574 

-0,146 

2,0 

+  0,107 

-0,043 

+0,418 

-0,274 

2,1 

+  0,087 

-0,057 

+0,340 

-0,363 

2,2 

+0,070 

-0,081 

+0,273 

-0,516 

3,0 

+0,043 

— 

+0,168 

— 

3,1 

+0,034 

-0,061 

+0,131 

-0,388 

3,2 

+0,015 

-0,060 

+0,058 

-0,381 
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Note  that  the  frequency  w  changes  somewhat  with  the  position  of  the 
measurement.  This  is  due  to  mass  loading  by  the  accelerometer.  The 
difference  is  sufficiently  small  to  be  insignificant,  and  is  assumed  to 
have  only  second-order  effects  on  the  mode  shapes,  in  view  of  the  wide 
separation  of  the  resonant  frequencies.  Finally,  note  that  is  a 
vector,  assumed  to  be  normal  to  the  blade  middle  surface.  The  direction 
cosines  relative  to  the  major  axis  of  the  root  were  measured  to  be: 


Station 

Inclination  a0 

cos  a 

0 

44° 

0.724 

1 

36 

0.808 

2 

31 

0.861 

3 

20 

0-938 

4 

13 

0.972 

4.  CALCULATION  OF  MODAL  INERTIA  AND  STIFFNESS  FROM  MODE  SHAPES 

The  thickness  of  the  blade  was  measured  at  several  stations  and 
across  the  chord.  Figure  92  shows  the  distribution  of  thickness.  The 
values  of  m^  and  fy  were  determined  from  Equations  45  and  46  using  the 
modal  data,  the  direction  cosines,  and  the  density  of  the  blade  material, 
assumed  to  be  0.27  lb/in3.  On  this  basis: 

m,  =  ff  ph  dxdy  -  0.235  kg 
xy 


0.235  x  1  x  1 

- - - - =  7.25 

0.0324 

ii^  =  °yy^  =  0.0324  kg 

.'.I  I  -  L^'iij  -  1.296  x  10  ’  kgnf 

i 

=  iryy  --  0.0324  {?  x  1 20 ) 2  =  18419  N/m 
K,  L?k'i  =  737  Nm/m 
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Figure  92.  Thickness  Distribution  of  Blade 
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Note  that  the  value  of  K-|  is  measured  in  N/m  because  X-j  is  the 
independent  variable,  not  X^/L,  where  L  is  the  blade  length. 

5.  MEASUREMENT  OF  RESPONSE  WITH  SLIP 

The  tests  with  slip  were  conducted  using  a  magnetic  transducer  for 
excitation,  with  a  pre-selected  gap,  (0.04  in.,  1mm)  at  several  different 
force  levels  ranging  from  0.0064  Newtons  to  0.125  Newtons,  the  maximum 
that  could  be  generated  with  the  transducer  available.  The  net  normal 
force  provided  by  the  two  springs  to  the  wires  at  the  root  of  the  blade 
was  varied  from  10  kg  to  25  kg.  As  the  excitation  frequency  was  increased 
slowly,  at  low  force  levels,  the  classical  shape  of  the  response  curve 
(|X^/S|  versus  frequency)  was  reproduced,  but  as  the  exciting  force  was 
increased  the  behavior  became  increasingly  nonlinear.  The  results 
obtained  for  tests  with  air  and  water  as  the  working  fluid  around  the 
blade  root  mating  surfaces  are  shown  in  Figures  93  to  1C9  for  the  first 
three  modes.  Test  results  are  summarized  in  Tables  6  to  10. 

Several  significant  facts  may  be  noted  from  these  figures.  One  is 
that  as  S  is  decreased,  a  point  is  reached  at  which  slip  never  occurs 
and  the  behavior  is  then  linear.  In  that  case,  the  only  significant 
remaining  source  of  damping  is  hysteresis  of  the  blade  material  itself. 

The  second  is  the  fact  that  as  water  is  substituted  for  air  as  the 
working  fluid  at  the  root,  slip  occurs  far  more  readily.  This  is 
accounted  for  by  the  fact  that  coefficient  of  friction  is  lower  in  this 
case.  Figures  101  and  102  show  some  of  the  acceleration  waveforms 
observed  on  the  oscilloscope  screen  for  various  conditions.  The  upper 
curves  represent  the  blade  acceleration,  and  the  lower  curves  the  force, 
all  at  the  point  of  maximum  amplitude.  In  air,  one  can  see  significant 
deviations  from  a  true  sinusoidal  shape.  In  water,  the  lower  threshold 
of  slip  seems  to  allow  high  order  harmonics  to  be  present  in  the  signal. 
The  reason  for  this  is  not  known. 
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Figure  99.  Measured  Response  in  Mode  2 
(N  =  10  kg.  Water) 


MO  *00  920  940  960  980 

FREQUENCY  Mi 


Figure  100.  Measured  Response  in  Mode  3 
(N  =  10  kg.  Air) 


TABLE  6  (CONCLUDED) 
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TABLE  7 


TEST  RESULTS  FOR  BLAOE  ROOT  IN  AIR  N  =  25  kg  (55.0  lb), 
CALIBRATION:  7.5  db  =  1  g  SPRING  TYPE  1,  LENGTH  =  24.1 
(0.95  INCHES) 
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TEST  RESULTS  FOR  BLADE  ROOT  IN  WA 
CALIBRATION:  7.5  db  =  1  g  SPRING 


TABLE  9  (CONCLUDED) 


139  -10  2.68 

149  -15  1.31 

170  -20  5.67E-! 

203  -24  2.51 


TABLE  10  (CONCLUDED) 
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Finally,  the  experiments  are  compared  with  analysis  in  Figures  103 
to  107. 

6.  DIRECT  MEASUREMENT  OF  MODAL  MASS  AND  STIFFNESS 

Section  2.3  describes  the  linear  modal  analysis  of  blade  response 
without  slip  (fully  clamped  root).  From  this  analysis  the  observed 
resonant  peak  acceleration  at  the  point  0,0  of  the  blade,  when  excited  at 
the  same  point  by  a  force  S,  is: 


c,  2 

i „  i  SL  cosa 

'  *1 '  T  2 

Ilnl“l 

2  2 

I  =  SL  cosa  _  SL  cosa  (127) 

I  X-|  |  n-j  |  X-j  [ 

From  Figures  93  to  97  for  the  fundamental  mode,  S  and  | (  at  low  force 
levels,  where  slip  does  not  occur,  are  known,  is  the  observed  modal 
damping  and  u>,  is  the  observed  resonant  frequency.  Table  11  gives  the 
collected  results.  It  is  seen  that  I -pi. 04  x  10  kgm  and  K^(=I^w^  )  = 

635  Nm/m.  These  values  are  very  close  to  those  determined  from  the 
measured  mode  shapes. 

7.  MEASURED  MODAL  DAMPING 

Further  understanding  of  the  test  results  is  afforded  by  examining  the 

f 

apparent  modal  damping  =  Af/f-j ,  where  f-j  is  the  resonance  frequency  and 
Af  is  the  frequency  3  db  below  peak  amplitude  (3.01  db  actually).  Clearly, 
since  the  system  behavior  is  highly  nonlinear,  this  linear  system  repre¬ 
sentation  is  not  particularly  appropriate  except  as  a  crude  numerical 

t 

indicator  of  the  effect  of  slip.  Even  so,  a  graph  of  versus  I/B2  is 
most  informative,  as  Figure  108  shows.  It  appears  that  the  graph  of  n-j 
versus  S/pN  is  nearly  a  straight  line,  for  each  fluid,  terminating  with 
a  discontinuity  at  n-j  =  0.01  (the  linear  damping)  and  finite  values  of 
S/uN.  The  agreement  between  analysis  and  experiment  is  good. 


Ill 
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FREQUENCY  Hz 


Figure  107.  Measured  and  Predicted  Response  for  Mode  2 
(S  =  0.084  Newtons,  N  =  10  kg.  Air) 
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Figure  108.  Measured  and  Predicted  Damping  Versus  1/6 
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Finally,  for  very  low  load  levels  the  residual  damping  is  linear  and 
apparently  hysteretic  in  nature,  and  its  source  is  possibly  the  blade 
material  itself.  For  the  first  mode  n-j  =  0.01  at  about  120  Hz.  Further 
tests  on  the  second  and  third  modes  show  that  =  0.006  and  113  =  0.0015 
for  the  second  and  third  modes,  at  420  Hz  and  920  Hz,  respectively. 

This  variation  of  nn  with  frequency  is  again  fairly  consistent  with 
published  data  for  such  damping. 

8.  CALIBRATION  OF  TEST  SYSTEM 

In  order  to  be  able  to  conduct  quantitative  measurements  of  the 
response  of  the  blade  to  a  known  harmonic  force,  it  is  necessary  to 
calibrate  the  measurement  and  excitation  systems.  For  example,  the 
voltage  across  the  fixed  resistor  in  the  exciter/oscillator  system,  at 
the  load  terminals,  must  be  monitored  and  related  in  some  way  to  the 
force  applied  to  the  blade.  Similarily,  the  signal  from  the  acceler¬ 
ometer,  after  passing  through  the  line  driver,  the  charge  amplifier  and 
the  various  cables  is  read  at  the  voltmeter,  and  the  relationship 
between  this  voltage  and  the  actual  acceleration  must  be  determined  for 
the  particular  system  and  components  used.  This  is  the  aim  of  the 
calibration  procedures,  which  must  be  conducted  prior  to  testing  and  from 
time  to  time  thereafter. 

The  accelerometer- Zero  Drive-Charge  Ampl i fier-Vol tmeter  measurement 
chain  is  calibrated  by  introducing  a  known  harmonic  acceleration  at  the 
accelerometer  and  simply  observing  the  R.M.S.  (root  mean  square)  voltage 
registered  at  the  voltmeter.  The  known  acceleration  is  generated  by  a 
small  electrodynamic  shaker,  having  a  steel  ball  embedded  in  a  hollow 
cavity  in  the  shaker  table,  as  illustrated  in  Figure  109  (step  1).  At  an 
acceleration  of  just  over  1  g  at  the  table,  the  steel  ball  will  separate 
from  the  floor  of  the  cavity  and  impacts  occur.  These  are  easily  noted 
as  "hash"  on  the  otherwise  smooth  sinusoidal  trace  on  the  oscilloscope 
screen.  The  acceleration  at  the  shaker  table  is  adjusted  until  the 
"hash"  just  appears.  Hence,  the  accelerometer  is  calibrated  to  1  g.  For 
the  Endevco  22  accelerometer,  MB  Zero  Drive  and  Amplifier,  and  HP  3400A 
voltmeter,  this  calibration  was: 

1  g  :  1.91  volts  (8  db) 
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Fiqure  109.  Calibration  Setups 


For  the  Endevco  Type  2226  accelerometer  used  in  the  calibration  of  the 
transducer,  the  appropriate  calibration  was: 

1  g  s  7.60  volts  (17.6  db) 

In  order  to  calibrate  the  transducer,  it  was  necessary  to  measure 
the  force  generated  on  a  block  of  iron  placed  a  distance  h  (the  gap) 
from  the  active  end  of  the  transducer  and  having  a  certain  oscillating 
voltage  V  at  the  output  terminals,  corresponding  to  a  definite  current  in 
the  circuit.  A  three  step  process  is  necessary,  since  the  force  cannot  be 
measured  directly.  Figure  109  shows  the  three  steps  in  this  sequence. 

The  first  step  is  to  calibrate  the  accelerometer,  as  discussed 
already.  When  this  step  has  been  accomplished,  one  has  a  numerical 
relationship  between  the  acceleration  seen  by  the  accelerometer  and  the 
voltage  registered  by  the  voltmeter,  in  volts  or  decibels  with  0.707 
volts  T  0  db,  because  this  is  an  R.M.S.  voltmeter. 
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The  second  step  was  to  calibrate  the  impedance  head,  so  that  it  in 
turn  can  be  used  to  measure  force.  This  was  accomplished  by  applying  a 
known  acceleration,  using  the  results  of  step  1 ,  to  a  shaker  table 
supporting  the  impedance  head.  The  force  experienced  by  the  crystal  of 
the  impedance  head  is  equal  to  (M  +  AM)X,  the  d'Alembert  force  corre¬ 
sponding  to  the  acceleration  X  acting  on  the  mass  MQ+AM  above  the  crystal. 
Mq  is  the  built-in  mass  of  the  impedance  head  itself  and  AM  is  the  added 
mass  of  the  soft  iron  piece.  Since  MQ  may  not  be  known,  one  can  vary 
the  mass  AM  and  measure  the  output  voltage  of  the  force  gage, 
registered  at  the  voltmeter  for  each  case.  With  X  known  from  the 
corresponding  accelerometer  reading,  we  can  then  plot  a  graph  of  V^/X 
against  Am  for  the  particular  force  gage  used  (Wilcoxon  Z  602).  The 
intersection  of  the  line  defined  by  the  measured  points  with  the  negative 
AM  axis  corresponds  to  the  mass  Mq  of  the  impedance  head  above  the  crystal. 
With  Mq  then  known,  any  reading  of  (MQ  +  M)X  gives  the  force  corresponding 
to  the  acceleration  X  and  the  voltage  output  from  the  force  gage,  via 
the  charge  amplifier,  the  zero  drive  and  the  cables.  It  was  found  that 
Mq  =  30  grams.  Hence  one  has  the  calibration;  for  example,  for  Vs  =  13.5 
volts.  Am  =  11.5  grams,  the  acceleration  X  was  0.168  g's.  Therefore: 

13.5  volts  3  (11,5-+  30)  x  0  168  x  9>81  Newtons 

1000 

or 

1  volt  =  5.08  x  10  ^  Newtons  (128) 

This  result  is  dimensioned  properly  by  converting  Mq  and  AM  to  kilograms 
and  the  acceleration  (0.168  g's)  to  meters  per  second. 

The  third  step  is  to  calibrate  the  transducer  itself,  using  the  force 
gage.  The  force  gage  was  placed  on  a  rigid  block  and  the  transducer  was 
then  brought  close  to  the  iron  block,  with  separation  h  as  shown  in 
Figure  109,  step  3.  The  voltage  V,  corresponding  to  a  particular  level 
of  the  oscillator  current  in  the  transducer  coils,  is  set  at  a  particular 
value,  such  as  30  volts.  The  gap  is  set  at  a  particular  value,  in  this 
case  0.043  inches  (1.092  mm),  by  means  of  a  feeler  gage  or  a  metal  sheet 
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of  this  thickness.  The  force  gage  was  then  connected,  again  through  the 
same  zero  drive,  charge  amplifier,  and  cables  to  the  voltmeter.  From 
the  force  gage  calibration,  given  in  Equation  128,  the  force  was  then 
calculated  directly  from  the  observed  voltage.  The  test  data  supports 
the  assumption  that  S  o  V,  but  only  one  gap  setting  was  tested  so  that 
the  other  assumption,  namely  S  a  1/h,  was  not  verified.  However,  the  same 
gap  was  used  in  all  tests,  so  this  omission  was  not  too  important.  The 
variation  of  Sh/V  with  frequency  is  shown  in  Figure  110.  One  can 
therefore  estimate  S  for  any  given  frequency  and  voltage  for  the 
particular  gap  setting. 

The  effective  centrifugal  load  on  the  blade  was  provided  by  two 
springs.  This  force  was  determined  by  measuring  the  change  in  length  of 
the  springs  under  load,  and  this  in  turn  was  calibrated  by  applying 
known  static  loads  to  each  spring  in  turn  and  measuring  the  corresponding 
length. 

9.  TESTS  ON  BLADE  WITH  SLIP  AT  A  PLATFORM 

With  the  specific  blade  used  in  the  earlier  experiments  on  blade  root 
damping,  the  use  of  a  gravity  loaded  rod  in  a  guide  allowed  a  controlled 
frictional  load  to  be  applied  at  selected  points,  as  illustrated  in 
Figures  111  and  112.  The  same  magnetic  transducer  and  pickup  accelerometer 
were  used  as  in  the  previous  tests,  and  the  gap  between  transducer  and 
blade  was  set  at  0.040  inches  (1.02  mm)  as  before.  Hence,  all  calibration 
factors  were  unchanged.  The  distance  X  was  2.25  inches  (57.2  mm),  and  at 
this  point  <})^  ( X/L )  =  0.20,  as  shown  in  Figure  90.  The  blade  was  excited 
at  various  force  levels  and  for  various  normal  loads  as  before.  The 
results  are  summarized  in  Table  12.  Figure  113  shows  typical  plots  of 
response  versus  frequency  for  the  undamped  blade  and  for  a  normal  load  of 
1.70  kilograms.  For  purposes  of  comparison  with  the  analysis  of  Section  3, 
we  note  that  the  mode  shape  in  the  fundamental  mode  would  not  change  too 
much,  so  that  M-j  -  6.036  kilogram  as  before.  For  the  blade  without  the 
normal  load,  the  stiffness  Kg  is  the  same  as  before,  i.e.,  Kg  =  22,500  N/m. 
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The  change  in  resonant  frequency  from  about  125  Hz  (unloaded)  to  148  Hz 
loaded,  leads  to  =  30,600  f'/m.  The  value  =0.2  corresponds  to 
the  modal  amplitude  ^ (5. 72/20)  -  0.20.  It  is  seen  that  analysis  and 
experiment  agree  quite  well.  A  friction  coefficient  p  =  0.3  was 
determined  for  slip  between  the  rod  and  the  blade  leading  edge. 


IO  100  1000  IOOOO 

FREQUENCY  -  Hi 


Figure  110.  Calibration  Curve  for  Transducer 
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Figure  112.  Test  Fixture 
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TABLE  12 


TESTS 

ON  BLADE 

WITH  SLIP 

AT  "PLATFORM" 

MODE  1, 

X  =  2.25 

in  (57.15 

mm),  R-j  =  0.20 

Input 

N 

Force  S 

al 

Freq 

Response 

X1 

(volts) 

(kg) 

Newtons 

s 

(H7.) 

(volts) 

(metres) 

40 

0 

0.079 

0 

127 

3.75 

7.70E-5 

126 

2.5 

5.21 

129 

2.5 

4.98 

130 

2.0 

3.92 

133 

1.0 

1.87 

139 

0.50 

8.57x15' 

159 

0.20 

2.6  2 

194 

0.10 

8.80x10' 

171 

0.15 

l.DxlO' 

152 

0.25 

3.59 

127 

4.0 

8.22E-5 

125 

2.0 

4.24 

123 

1.0 

2.19 

118 

0.5 

1.19 

107 

0.2 

5.79E-6 

95 

0.1 

3.67 

74 

0.0  5 

3.03 

48 

0.023 

3.31 

20 

0.017 

40 

0.02 

40 

1.70 

0.079 

11.2 

147 

1.50 

2.3DE-5 

143 

1.05 

1.70 

152 

1.05 

1.51 

153 

1.00 

1.42 

162 

0.50 

6.31E-6 

185 

0.20 

1.94 

222 

0.10 

6.72F.-7 

40 

3.52 

0.079 

23.1  149 

3.20 

4.78E-5 

146 

2.40 

3.73 

151 

2.40 

3.49 

152 

2.00 

2.87 

155 

1.00 

1.38 

164 

0.50 

6 . 16E-6 

191 

0.20 

1.82 

226 

0.10 

6.49E-7 

147 

4.00 

6.13E-5 

146 

3.00 

4.66 

146 

2.00 

3.11 

143 

1.00 

1.62 

138 

0.50 

8.70E-6 

123 

0.20 

4.38 

103 

0.10 

3.12 

82 

0.05 

2.46 

43 

0.02 

3.58 
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TABLE  12  (CONCLUDED) 


Input 

(volts) 

N 

kg 

Force  S 

N 

"l 

S 

Freq 

(HZ) 

Response 

(volts) 

X1 

(m) 

40 

2.61 

0.079 

17.1 

148 

2.50 

3.78E-5 

145 

1.75 

2.76 

151 

1.75 

2.54 

156 

1.00 

1.36 

164 

0.50 

6.16E-6 

194 

0.20 

1.76 

226 

0.10 

6.49E-7 

147 

2.60 

3.99E-5 

145 

1.50 

2.36E-5 

142 

1.00 

1.64 

137 

0.50 

8.83E-6 

122 

0.20 

4.45 

100 

0.10 

3.31 

79 

0.05 

2.65 

63 

0.03 

2.50 

143 

1.00 

1.62E-5 

137 

0.50 

8.83E-6 

123 

0.20 

4.38 

104 

0.10 

3.06 

83 

0.05 

2.40 

65 

0.03 

2.35 

40 

1.35 

0.079 

8.84 

145 

2.10 

3.31E-5 

143 

1.50 

2.43 

150 

1.50 

2.21 

154 

1,00 

1.40 

163 

0.50 

6.24E-6 

189 

0.20 

1.86 

225 

0.10 

6.54E-7 

146 

2.20 

3.42 

142 

1.00 

1.64 

136 

0.50 

8.96E-6 

121 

0.20 

4.53 

102 

0.10 

3.18 

81 

0.05 

2.52 

62 

0.03 

2.59 

40 

1.70  0.079 

11.2  145 

1.50 

2.35E-5 

142 

1.00 

1.64 

137 

0.50 

8.83E-6 

122 

0.20 

4.45 

100 

0.10 

3.31 

79 

0.05 

2.65 

63 

0.03 

2.50 

0.0405 
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Figure  113.  Response  of  Blade  with  Slip  at  Platform 
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SECTION  VI 
CONCLUSIONS 

The  analytical  and  experimental  investigations  described  in  this 
report  have  shown  that  the  apparent  modal  damping  in  a  simple  dovetail 
blade  can  be  predicted  on  the  basis  of  assumed  gross  slip  at  the  blade 
root-disc  interface,  provided  that  the  exciting  force  exceeds  a  certain 
threshold  amplitude.  For  driving  forces  below  this  level,  the  damping 
arises  from  other  sources  such  as  hysteresis  in  the  blade  material  or 
non-slip  interaction  at  the  blade  root-disc  interface.  The  investigation 
has  been  restricted  to  a  relatively  simple  problem,  and  much  more  effort  is 
needed  if  one  is  to  design  and  build  blades  which  are  optimized  with 
respect  to  slip  damping  at  high  rotation  speeds.  However,  the  basic 
analytical  techniques  and  physical  modelization  remain  the  same. 

It  has  also  been  shown  that  a  specific  configuration  of  compressor 
or  turbine  blade  root  geometry,  in  which  part  of  the  root  contacts  the 
disk  only  at  high  rotational  speeds,  can  provide  high  levels  of  slip 
damping,  provided  that  the  relevant  stiffnesses  are  properly  selected. 

The  analysis  can  provide  the  basis  for  preliminary  design  investigations, 
but  laboratory  and  spin  pit  testing  will  be  necessary  to  establish  the 
accuracy  of  the  approach  and  to  develop  the  specific  configurations  most 
appropriate  for  practical  application.  The  changes  in  blade  geometry 
needed  to  optimize  this  type  of  damping  do  not  represent  very  large 
departures  from  current  practice,  nor  need  they  represent  any  weight 
increases  over  current  blades. 

Apart  from  experimental  and  spin  pit  testing,  a  logical  follow  on 
effort  would  investigate  the  effects  of  finite  disk  compliance,  and  hence 
of  multiple  blade  systems,  on  the  dynamic  behavior.  However,  one  would 
expect  high  levels  of  slip  damping  to  still  be  attainable.  It  is  hoped 
that  this  report  will  stimulate  or  encourage  such  investigations,  since 
the  need  for  high  damping  in  rotating  blades  is  becoming  ever  more 
urgent. 
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